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Abstract 

The  source  of  solid  rocket  motor  vibration  is  investigated  by  treating 
the  problem  of  pressure  fluctuations  inside  the  motor  cavity.  Theoretical 
aeroacoustics  is  applied  to  the  cavity  gases  and  imbedded  burning  metal 
agglomerates.  Several  critical  experiments  are  performed  to  provide  nu- 
merical input  to  the  theory.  Consideration  of  turbulence,  combustion  and 
entropy  noise  yields  the  conclusion  that  only  ore  cause  is  dominant  for 
the  fluctuations  in  chamber  pressure  - that  of  interaction  of  turbulence 
with  the  exhaust  nozzle.  Typically,  an  0.4%  rms  pressure  fluctuation  can 
be  accounted  for  by  this  mechanism.  Spectral  distributions  of  the  noise 
are  presented  and  a comparison  is  made  of  the  theory  and  an  actual  motor 
firing. 
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Nomenclature 


A feedback  pacemeter 

a radius 

a^  isentropic  admittance  coefficient 

entropy  admittance  coefficient 
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Cp  specific  heat  at  constant  pressure 

it 
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d diameter 
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drag  in  i direction 
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1 length 
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m mass  addition  rate  per  unit  volume  or  particle  mass  or  mass 

flow  rate  per  unit  area 

57;  combustion  noise  source  function 

M Mach  number 

n bum  rate  pressure  exponent  or  number  of  droplets  per  unit  volume 

p pressure 

P perimeter 
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radial  coordinate  or  bum  rate 

Reynolds  number 

entropy 

area,  correlation  area 
Strouhal  number 
time 

sample  time  for  Fourier  transform 
temperature 


. th 


u velocity  in  i direction,  velocity  in  axial  direction 
th 

velocity  in  i direction 
1/n 

chamber  volume,  correlation  volume 
coordinate  in  i*"^  direction,  axial  coordinate 
i^^  coordinate  direction,  axial  coordinate 
mass  fraction  of  droplets 

droplet  drag  function  defined  in  Eq,  (11)  of  Appendix 

i 

mean  square  value 

thermal  diffusivity  or  coefficients  in  boundary  condition 
or  droplet  parameter  of  Eq,  (7)  of  the  Appendix 

isentropic  specific  wall  admittance,  boundary  condition 
parameters  or  droplet  parameter  of  Eq,  (8)  of  the  Appendix 

ratio  of  specific  boats 

source  function  of  Eq,  (12)  of  the  Appendix 
or  Dirac  delta  function 

coordinate  in  i^^  direction 
angular  coordinate 


complex  wavenumber 
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\ 

IJ- 

p 

a 


defined  in  Eq,  (12) 
gas  viscosity 
feedback  parameter 
coordinate  in  i*"^  direction 
density 

dimensionless  entropy  fluctuation,  s^/c^ 


m circular  frequency 

Q frequency  parameter  in  feedback  law 


Subscripts  except  where  defined  above 
I ,1.  ,I  ,11  ,11.  , III  ,IV  ,IV.  source  term  identifiers 

H D C d D AD 

d condensed  phase  material 

e nozzle  entrance  plane  or  effective  value 

h head  end 

inf  influence 

i,J  cartesian  coordinate  directions 

max  iaaxir.ium 

min  minimum 

nit  nitrogen  in  pressurization  vessel 

p pressure  perturbation  or  particles 

prop  immediately  above  propellant  surface 

r radial  component 

rms  root  mean  square 

5 condensed  phase 

w propellant  surface 

oo  plane  wave  mode 

6 pertaining  to  6 
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u)  Fourier  transform, 

Superscripts  except  where  defined  above 
' perturbation  quantity 

~ coefficient  in  boundary  conditions 

mean  value 
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I,  Introduction 


Hie  propulsion  community  has  traditionally  resigned  itself  to 
acceptance  of  a background  noise  level  in  rocket  motor  pressure  of 
the  order  of  1%  of  the  mean  chamber  pressure.  There  has  been  little 
concern  over  the  cause  of  such  fluctuations,  or  their  reduction,  ex- 
cept in  :ases  where  a)  they  have  substantially  exceeded  the  1%  level 
or  b)  taey  have  become  phase  coherent  instabilities  in  the  motor  or 
engine.  These  chamber  pressure  fluctuations  are,  of  course,  a source 
of  motor  vibration.  There  are  reasons  why  it  would  be  desireable  to 
reduce  this  vibration  level.  Consequently,  this  program  was  und“r- 
taken  to  understand  the  origin  of  the  chamber  pressure  fluctuations 
and  to  understand  the  scaling  rules  in  order  to  be  able  to  design 
low  vibration  into  a motor. 

This  report  is  concerned  with  solid  rocket  motors.  Potential 
sources  of  noise  in  the  chamber  pressure  are  relatively  easy  to 
identify,  although  quantification  is  much  more  difficult.  First, 
in  the  case  of  composite  solid  propellants  the  propellant  is  heter- 
ogeneous in  structure,  Thio  forces  an  unsteady  combustion  process, 
and,  hence  fluctuations  in  chamber  pressure.  Secondly,  the  flow  is 
turbulent  within  the  motor  cavity.  There  are  pressure  fluctuations 
due  to  the  turbulence  itself,  there  may  be  a turbulence-combustion 
interaction  called  direct  combustion  noise, and  there  may  be  an 
interaction  of  the  turbulence  with  the  nozzle  flow  called  nozzle- 
vorticity  interaction  noise, Thirdly,  there  may  be  hot  spots  in 
the  flow  field  due  to  metal  burning  in  the  chamber  volume;  these 
hot  spots  (or  cold  spots)  encountering  the  nozzle  will  give  rise  to 
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(3) 

noise  which  has  been  called  "entropy  noise",'  ' 

To  complicate  the  problem,  there  will  be  feedback  between  pressure 
fluctuations  and  the  combustion  process,  and  there  will  be  damping  of 
acoustic  energy  due  to  the  presence  of  condensed  phase  material  in  the 
chamber  gas  flow.  Moreover,  there  will  be  damping  due  to  the  exhaust  of 
acoustic  energy  through  the  nozzle.  The  chamber  pressure  fluctuation 
which  is  ultimately  achieved  will  depend  upon  a balance  of  the  acoustic 
losses  and  the  strength  of  the  sources  in  the  presence  of  feedback  with 
tne  combustion  process.  This  report  attempts  to  calculate  the  stationary 
random  oscillation  which  will  be  achieved  for  a given  motor  design.  The 
treatment  is  necessarily  approximate,  Wha'  is  primarily  sought  is  a)  the 
proper  order  of  magnitude  of  the  oscillation,  b)  a spectral  shape  esti- 
mate and  c)  scaling  rules  for  the  pressure  fluctuation  magnitude. 

For  reasons  to  be  delineated  below  the  theory  will  be  based  on  the 
Lighthill  approach to  aeroacoustics,^^^  VIhere  data  are  required  to 
feed  into  the  theory,  appropriate  data  from  the  literature  will  be  used, 
if  it  exists.  In  some  areas,  however,  exploratory  experiments  are  re- 
quired to  generate  the  data;  these  experiments  are  carried  out  and  re- 
ported on  in  this  report. 
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II.  Results  of  the  Theory 


The  details  of  the  theoretical  formulation  and  solution  are  con- 
tained in  the  Appendix.  What  is  sought  is  an  expression  for  the  pres- 
sure fluctuation,  and  its  spectral  content.  Accounted  for  in  an  approxi- 
mate manner  are  particulate  matter  suspended  in  the  gases,  the  effect 
of  a mean  flow  and  leedback  between  the  pressure  fluctuations  and  com- 
bustion process.  These  are  effects  which  are  important  in  the  chamber 
acoustics.  The  sources  of  the  pressure  fluctuations  are  the  turbulence, 
turbulence-combustion  interaction,  and  turbulence-nozzle  interaction. 

All  of  these  have  their  representations  in  the  theory. 

The  theory  specifically  considers  a stationary  random  oscillation. 

At  sufficiently  high  values  of  the  propellant  feedback  response  function 
the  motor  (and  the  computation)  will  indicate  an  instability.  At  this 
point  the  analysis  will  fail.  The  motor  must  be  stable  for  this  calcula- 
tion to  be  valid.  Even  in  the  case  of  a stable  rocket  motor,  there  is 
a background  noise  level  of  pressure  fluctuations.  The  calculation  of 
the  background  level  is  the  purpose  of  this  program. 

Major  limitations  of  the  theory  are  a)  frequencies  are  restricted 
to  those  lelow  cut-on  of  the  first  transverse  mode  of  oscillation  of  the 
chamber  gases  and  b)  the  cross-section  area  of  the  cavity  must  be  reason- 
ably invariant  with  axial  position.  The  effects  of  these  restrictions 
are  that  a)  only  axial  vibration  of  the  motor  is  under  consideration 
and  b)  there  may  be  difficulties  in  interpretation  of  data  from  motors 
with  a complex  grain  configuration.  The  theory,  therefore,  only  considers 
nearly  plane  wave  motion  in  the  motor  cavity  so  that  the  pressure  fluctua- 
tion, p%  is  a function  of  axial  length,  x,  alote.  Vibration  of  a motor 
in  the  axial  direction  is  set  up  by  the  difference  in  pressure  acting  on 
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the  head  end  as  opposed  to  the  tail  end,  idiere  the  nozzle  is  discharging 
the  gases  and  particles.  Roughly,  the  thrust  fluctuation  is  F^as  Aj.  so 
tha  pressure  fluctuation  is  the  primary  quantity  of  interest. 

The  mean  square  pressure  fluctuation,  < p^^>  is  related  to  its 
spectral  distribution,  Gp(f)  by^^^ 


< p'2  > = J G (f)df 


(1) 


where  Gp(f)  may  be  estimated  by  a Fourier  transform  operation  on  p^  by 

P * P 

G (f)  = 2 ^ 

P ’^o 


p'(t)e 


•2  nif 


dt 


(2) 


Consequently,  the  solution  to  the  aeroacoustic  problem  is  given  in 
terms  of  the  Fourier  transform  of  the  pressure  fluctuation  as  the 
quantity  of  interest. 

As  shown  in  the  appendix  the  solution  for  G^  at  the  boundary  of 
the  chamber  (say,  the  head  end)  is  given  by 


G 

P 


= G + 

Pi 

a 


+ G 


+ G 


+ G 


+ G 


"Ill 


(3) 


where  each  term  in  the  sum  is  due  to  a specific  noise  source.  In  Eq. 

(3)  an  assumption  is  that  there  is  no  correlation  between  the  different 
noise  sources.  This  assumption  is  open  to  some  question  but  becomes 
unimportant  in  the  end,  because  only  one  of  the  terms  of  Eq.  (3)  is  found 
to  be  important.  In  Eqs,  (3)  the  terms  designated  by  I have  to  do  with 
noise  caused  by  pure  turbulence,  I and  I,  are  related  to  the  turbulence 

d D 

distributed  through  the  cavity  volume  and  are  separated  into  two  terms 
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purely  for  analytical  convenience.  Term  is  caused  by  the  interaction 

of  the  turbulence  and  the  nrzzle  flow  process.  Term  II  deals  with  noise 

3i 

generated  by  combustion  near  the  surface  of  a heterogeneous  propellant 
and  term  IIj^  is  combustion  noise  of  particulate  metal  imbedded  in  the 
gas  stream.  Term  III  is  noise  from  an  interaction  of  hot  spots  with 
the  nozzle  flow  process. 

Written  out  in  full  the  individual  terms  of  Eq.  (3)  become 


Each  term  of  Eqs,  (4)  contributes  to  the  square  of  the  pressure  fluc- 
tuation, so  a square  of  a source  term  occurs  in  each  of  the  above. 

In  each  case  either  the  volume  or  surface  integral  contains  the  number 

g as  a factor  in  the  integrand.  This  is  the  Green's  function  for  the 
U) 

duct  acoustics  problem  and  tells  how  a source  gets  turned  into  pressure, 
depending  upon  the  location  and  the  strength  of  the  source.  The  Green's 
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function  and  its  behavior  are  described  in  the  Appendix,  To  calculate 
it  requires  the  flow,  particulate  and  nozzle  details  and  also  the 
behavior  of  the  feedback  response  law  for  the  propellant.  At  low 
frequencies  an  analytical  expression  for  the  Green's  function  is 
worked  out  in  the  Appendix,  It  is  valid  for  frequencies  lower  than 
about  1/2  the  first  longitudinal  mode  frequency.  The  expression  for 
the  square  of  the  Green's  function  is 


Vu)  ^2  g2  _ ^^^2 


(5) 


This  is  useful  for  quick  order  of  magnitude  use  in  Eqs.  (4)  and  for 

rough  investigation  of  scaling  rules. 

In  each  of  the  terms  of  Eqs,  (4)  tliere  is  a source  term  to  be 

evaluated.  For  example,  in  G^  , there  is  required  knowledge  for  T^j^, 

^a 

which  is  a turbulence  fluctuation  term.  The  estimation  of  these 
terms  is  deferred  until  Section  IV, 

Although  it  was  not  guaranteid  at  the  outset,  the  theory  has 
recovered  every  physically  expected  noise  source,  VJhile  the  accuracy 
might  not  be  high,  because  the  simplest  aeroacoustics  formalism  was 
used  and  several  approximations  were  introduced,  the  scaling  rules 
and  order  of  magnitude  of  various  terms  should  be  correct.  Moreover, 
the  accuracy  with  which  the  input  variables  are  known  i;  low,  as  will 
be  seen,  because  of  the  lack  of  some  fundamental  data  on  turbulence 
in  solid  rocket  motor  cavities.  Consequently,  : t is  considered  that 
the  accuracy  of  the  theory  and  the  accuracy  of  the  input  data  are 
commensurate. 
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Ill,  Experimental  Results 


After  viewing  the  results  of  the  theory,  in  Eqs,  (4),  there  are  several 
experimental  inputs  required.  These  are  spatial  correlation  length  scales 
and  spectral  behavior  of  the  axial  velocity  fluctuations,  temperature  fluc- 
tuations and  mass  flow  per  unit  area  fluctuations  just  above  the  propellant 
surface.  None  of  these  quantities  have  been  investigated  in  the  literature 
in  the  detail  required  for  this  program.  However,  reasonable  magnitude  esti- 
mates may  be  made  of  the  axial  velocity  turbulence  characteristics  through 
use  of  Laufer's  pipe  flow  data^^^as  corrected  by  the  results  of  Ref,  8, 
which  investigated  the  effects  of  mass  addition  from  the  wall.  The  other 
data  is  generated  here  through  experiments  which, although  not  precision 
experiments,  will  yield  preliminary  magnitude  and  spectral  estimates  for 
the  unknown  quantities, 

A,  Combustion  Noise  from  Heterogeneous  Propellant 

The  object  of  this  measurement  was  to  determine  the  strength  of  the 

combustion  noise  source  if  the  propellant  were  radiating  sound  to  a free 

field.  This  information  can  be  translated  to  knowledge  of  in  Eqs,  (4), 

The  apparatus  shown  in  Fig,  III  A-1  was  used  to  determine  the  acoustic 

properties  of  the  burning  propellant.  This  tube  had  a length  of  64"  (162,5cm) 

and  an  inside  diameter  of  4"  (10,16  cm),  It  was  described  in  detail  in  Ref.  (9), 

The  tube  had  a maximum  operating  pressure  of  1000  psi, 

A typical  sample  is  shown  in  Fig,  III  A-2,  Tlie  strands  of  propellants 

were  cut  to  be  V'  x x 3" (6,35  mm  x 6,35  mm  x 76,2  mm),  TSiis  size  of  sample 

resulted  in  an  average  action  time  pressure  approximately  10%  higher  than  the 

initial  pressure  xdien  the  tube  was  pressurized  to  300  psi.  The  sides  of  the 

strands  used  in  this  investigation  were  coated  with  approximately  l/8"(3,2ram) 

of  Dow  Corning  silicon  rubber  coating  to  act  as  an  inhibitor  to  the  spread 

of  deflagration  down  the  side  of  the  end  burning  strands.  This  thick  coating 
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Figure  III  A-2,  Propellant  Strand  Prepared  for  Acoustic 
Testing, 
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was  necessary  to  prevent  metal  agglomerate  recirculation  ignition  of  the 
sample  along  the  side  of  the  sample  when  aluminized  propellants  were  tested. 
Once  the  inhibitor  on  the  side  of  the  sample  was  cured  the  sample  was  trim- 
med and  mounted  on  an  aluminm  flange  using  the  same  silicon  rubber.  The 
ignition  power  supply  required  approximately  1 3/4"  (4.5  cm)  of  .010"(.25mm) 
diameter  nichrome  wire.  This  wire  was  formed  into  a grid  and  bonded  to  the 
V'  X V'  (6»35  mm  x 6.35  ram)  surface  with  a borau* titanium  based  ignitor 
paste  (x  - 225  pyrotechnic  mixture).  Two  6"  electrical  leads  of  #28  teflon 
insulated  copper  hookup  wire  were  soldered  to  the  ignitior  wire  for  hookup 
to  the  electrical  feed-through  in  the  acoustic  tube.  These  wires  were  then 
tied  to  the  side  of  the  inhibited  sample  to  minimize  the  chance  of  ignition 
wire  damage  during  handling. 

Hie  prepared  sample  was  bolted  on  the  one  inch  (2,54  cm)  thick  closure 
flange  of  the  acoustic  tube,  A special  viscous  grease  was  used  between  the 
aluminum  sample  flange  and  the  stainless  steel  closure  flange  to  ensure  a 
transmission  path  for  ultrasonic  acoustic  emissions.  Both  of  these  flanges 
had  precision  ground  surfaces  at  the  mounting  points,  A Dunegan-Endevco 
Model  D 9201  ultrasonic  transducer  was  mounted  on  the  opposite  side  of  the 
closure  flange  to  measure  the  ultrasonic  acoustic  emission  of  the  burning 
propellant.  This  signal  was  used  to  determine  (he  bum  rate,  Ihe  prepared 
sample  is  located  at  the  right  hand  end  of  the  acoustic  emission  tube  shown 
in  Fig,  III  A-1, 

An  acoustic  hard  termination  with  a vent  hole  in  the  center  was  located 
7, 87"  (20  cm)  away  from  the  left  hand  end  closure  flange.  The  acoustical 
transducer,  a PBN  Model  376A  or  376LF,  piezoelectric  sensor  was  located 
immediately  in  front  of  the  tiard  termination,  A CEC  Model  20077  pressure 
cransducer  was  locaced  in  the  exhaust  line  to  record  the  pressure  rise  due 
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to  combustion  of  the  solid  propellant  in  the  acoustic  tube. 

The  three  signals  were  filtered  and  amplified  and  then  recorded  on  a 
magnetic  tape  recorder.  The  acoustic  data  and  tlie  pressure  data  were  re- 
corded using  the  FM  mode  for  more  accurate  signal  reproduction.  Recording 
at  60  in/sec  (152,4  cm/sec)  gave  an  upper  limit  of  30,000  Hz,  It  was  nec- 
essary to  record  the  ultrasonic  signal  using  the  direct  or  AM  mode  of 
signal  recording.  This  gave  an  upper  limit  of  300  kHz,  Tht  data  acquistion 
system  is  shown  schematically  in  Fig,  III  A-3, 

The  signals  were  then  played  back  at  reduced  rates  for  graphical  re- 
production and  Fourier  analysis.  Some  signal  conditioning  in  the  form  of 
filtering,  discriminating  and  counting  was  necessary  for  all  channels.  The 
data  reduction  system  is  shown  schematically  in  Fig,  III  A-4, 

Typical  results  for  an  aluminized  propellant,  MC-170,  are  shown  in  Fig, 
III  A-5,  The  pressure  rise  for  this  test  was  from  300  psi  to  366  psi,  llie 
average  pressure  was  333  psi.  The  acoustic  signal  detected  by  the  acoustic 

transducer  represented  an  average  sound  pressure  level  of  91  db  re  2 x 10*"^ 

2 

N/m  • The  burn  time  indicated  by  the  ultrasonic  transducer  agreed  with  the 
pressure  rise  and  the  acoustic  signal.  It  was  a uniform  burning  strand  with 
an  initial  burst  of  ultrasonic  noise  due  to  the  ignitor  paste  and  a final 
burst  due  to  the  non-uniform  burning  as  the  deflagration  reaches  the  mount- 
ing flange  and  the  mounting  material.  This  sample  was  3,5"  (8.9  cm)  long. 
The  burn  rate  for  this  propellant  was  ,39  cm/sec  (1  cm/sec)  at  300  psi. 

The  pressure  spectral  level  of  the  acoustic  signal  was  obtained  using 
the  Hewlett  Packard  5451A  Fourier  Analyzer.  The  preliminary  analysis  of 
the  acoustic  signal  shown  in  Fig,  III  A-5  is  shown  in  Fig,  III  A-6,  The 
range  of  this  analysis  was  from  4 Hz  to  1024  Hz  with  a resolution  of  4 Hz. 
This  was  obtained  by  using  a block  size  of  512,  a time  window  of  2 seconds 
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Figure  III  A-3.  Data  Acquisition  System  for  Combustion 
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and  a speed  redaction  of  8:1  on  the  magnetic  tape  recorder.  The  signals 
were  digitized  and  two  Hanning  windows  were  used  to  reduce  leakage.  All 
acoustic  signals  obtained  in  this  section  of  the  investigation  were  fil- 
tered with  a Krohn-Hite  electronic  filter,  model  3323R,  set  to  pass  all 
frequencies  above  5 Hz.  This  was  necessary  to  prevent  saturation  of  Neff 
amplifiers  due  to  large  scale,  low  frequency  oscillations  of  the  acoustic 
signal,  niese  oscillations  were  initiated  by  the  ignition  and  starting 
transient  of  the  combustion  process.  The  magnitude  and  duration  of  these 
oscillations  were  enhanced  by  the  BBN  power  supply.  Model  P-16,  and  lower 
high  pass  frequency  settings  of  the  Krohn-Hite  filter.  A 60  Hz  noise  prob- 
lem in  the  recorded  acoustic  signal  was  discovered  when  the  data  reduction 
for  this  series  of  tests  was  initiated.  This  was  traced  to  the  Krohn-Hite 
filter.  This  60  Hz  signal  was  within  the  specifications  of  the  filter  and 
only  became  a problem  when  it  was  amplified  by  the  Neff  amplifier  at  a 
gain  of  500  or  1000.  The  occurance  of  the  60  Hz  signal  and  its  harmonics 
at  comparable  levels  with  the  combustion  generated  noise  is  shown  in  Fig, 

III  A-6,  Hiere  was  sufficient  recorded  signal  to  obtain  background  pressure 
spectral  levels  for  all  propellants  at  300  psi,  before  sample  ignition. 

This  signal,  labeled  background,  was  subtracted  from  the  uncorrected  spectra, 
giving  the  corrected  pressure  spectral  level.  This  reduced  the  peaks  in  the 
pressure  spectra  at  60,  180  and  300  Hz, 

These  corrected  pressure  spectra  were  further  reduced  to  yield  the  strength 
of  the  combustion  noise  source  if  the  propellant  were  radiating  sound  to  a free 
field.  The  corrected  pressure  spectral  level  shown  in  Fig.  Ill  A-6  can  be  re- 
lated to  the  actual  pressure  fluctuation  in  the  following  manner 

FSL/^H.  ' 10  V lvou2 
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This  Gp  can  be  related  directly  to  the  calibration  of  the  acoustic  trans- 
ducer. The  tape  recorder  had  a voltage  reduction  of  4,47,  the  Neff  amplifer 
had  a gain  of  500  and  the  acoustic  calibration  of  the  transducer  for  a sig- 
nal 160  db  re  .0002  dy/cm^  gave  ,122  volt,  RMS,  The  actual  strength  of  the 
combustion  noise  source,  G^,  as  required  in  the  theory, is  related  to 
the  pressure  spectral  level  G^  as  follows. 


“ic 


prop 


nit 


Gp  (sin  \l  + P^cos  kl) 


The  quantity  p was  evaluated  at  the  first  peak  and  trough. 


(10) 


G (peak)/-  = 1/p^ 

P ^p (trough) 

Three  propellants  were  tested  initially.  They  trere  MC-170,  UTP-3001  and 
a GMDB  propellant.  The  CWDB  propellant  was  obtained  from  the  U,  S,  Army 
BMDATC  for  this  program.  The  composition  of  these  propellants  along  with 
that  of  two  other  propellants  MC-172  and  T-48,  is  given  in  Table  III  A-1, 

The  sample  bum  rate,  tube  average  pressure  and  the  average  sound  pres- 
sure level  for  the  comparison  of  the  two  BBN  transducers,  models  376A  and 
376liF(  the  low  frequency  transducer)  are  shown  in  Table  III  A-2,  The  actual 
strength  of  the  combustion  noise  source,  G^,  is  shown  in  Fig,  III  A-7  for 
the  six  tests.  The  pressure  spectral  levels  were  corrected  for  the  frequency 
range  4 Hz  through  the  third  trough  of  Fig,  III  A-6, 

The  results  for  MC-170  using  the  low  frequency  transducer  do  not  agree 
with  the  other  two  propellant  results.  Uie  sample  did  not  bum  linearly. 

This  is  shown  in  Fig.  Ill  A-8  \^ere  the  pressure  signal  does  not  agree  with 
the  ultrasonic  signal.  The  rapid  pressure  increase  has  driven  the  acoustic 
transducer  to  saturation.  This  transducer  appears  to  recover,  but  there  is 
no  acoustic  signal  for  the  rest  of  the  run.  This  behavior  was  noted  for  all 
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Table  III  A-1  Properties  of  Propellants  used  in  Tests 


Designation 

AP 

A1 

Binder 

MC-170 

200p,in 

26.8% 

5|j,m  20% 

HTPB  12% 

14p,in 

20.4% 

68% 

6(j,m 

20.8% 

UTP-3001 

190nm 

44.3% 

38;xm  16.1% 

PBAN  16.15% 

8p,ra 

23.2% 

67.5% 

MC-172 

200^in 

26.8% 

5p,in  20% 

HTPB  12% 

14jj,m 

20.4% 

68% 

(AFCAM) 

6y,in 

20.8% 

T-48 

48p,m 

70% 

HTPB  30% 

C31DB 

FZOi 

Double 

Base,  HMX  56% 
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Table  III  A— 2 Results  of  Comparison  of 
Propellants  and  Transducers 


Propellant 

Transducer 

(BBN) 

MC-170 

376A 

MC-170 

376LF 

UTP-3001 

376A 

UTP-3001 

376LF 

(MDB 

376A 

CMBD 

376LF 

Burn 

Rate 

(in/sec) 

Average  Pressure 
(psi) 

.39 

333 

.52 

353 

.39 

338 

.37 

338 

.175 

328 

.179 

336 

Sound  Pressure  Level 

db  re, 0002  dy/  2 
cm 

91 

105 

93 

90 

78.5 

77.5 
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Figure  III  A-7.  Strength  of  the  Combustion  Noise 


tests  using  the  new  transducer  (376LF)  in  this  series.  All  Fourier  analysis 
was  performed  on  the  portion  of  the  trace  before  any  transducer  saturation 
occurred.  The  first  point,  at  4 Hz,  shown  in  Fig.  Ill  A-7  is  in  error  since 
the  signal  was  passed  through  a high  pass  filter  set  at  5 Hz  before  it  was 
amplified  and  recorded.  Additional  tests  were  scheduled  to  obtain  pressure 
spectra  at  lower  frequencies  (i.e,  less  than  5 Hz). 

A series  of  eight  direct  noise  tests  were  completed  using  both  MC  170 
and  MC172.  The  results  presented  here  are  for  MC  172,  idilch  is  very  similar 
to  MC  170  (see  Table  III  A-1),  Three  of  the  tests  were  conducted  using  vari- 
ous amplifier  gains  with  the  low  frequency  transducer,  BBN  376LF.  No  useful 
data  were  obtained  due  to  tape  recorder,  amplifier,  and  transducer  satura- 
tion for  all  values  of  gain  used.  The  next  four  tests  were  conducted  using 
another  type  of  piezoelectric  transducer,  AVL  type  8QP500  ca.  It  was  thought 
that  this  transducer  and  associated  charge  amplifier  would  have  a better  low 
frequency  response  and  recover  from  saturation  conditions  faster.  This  was 
not  obvious  after  the  Fourier  analysis  of  the  resultant  tests. 

The  final  test  of  this  series  used  the  BBN,  376LF,  low  frequency  trans- 
ducer with  a modified  power  supply.  The  current  regulating  diode  was  replaced 
with  a factory  supplied  diode  to  allow  for  higher  transducer  currents,  Hiis 
increased  battery  drain, and  it  did  improve  the  low  frequency  response  and  re- 
covery, The  pressure  spectral  level  is  shown  in  Fig,  III  A-9,  Hie  signal  due 
to  combustion  was  at  least  35  db  higher  than  the  backgi  .nd  measurements  over 
the  1,6  to  30  Hz  range.  This  separation  of  the  low  frequency  spectra  was  visi- 
ble at  reduced  gain  with  a high  pass  filter  blocking  the  DC  to  1 Hz  signal. 

It  was  also  present  in  the  last  test  using  the  AVL  transducer  but  the  separa- 
tion v\ried  from  5 to  10  db  over  the  1,6  to  50  Hz  region.  The  nujor  result. 
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the  very  low  frequency  data  was  obtained,  was  that  there  was  no  sur- 
prise at  low  frequency.  The  noise  spectra  level  off  as  zero  frequency  is 
approached.  No  further  direct  noise  tests  were  conducted,  since  the  magni- 
tude analysis  showed  that  the  overall  level  of  noise  was  too  low  to  be  im- 
portant to  the  vibration  problem  for  the  propellants  tested. 

B.  Cross-flow  Effects  on  Combus tion  Noise 

The  object  of  the  experiment  was  to  measure  the  combustion  noise  from  a 
process  where  there  was  a combustion  and  flow  interaction.  This  effect  has 
been  explicity  excluded  from  the  theory  and  it  was  desired  to  know  the  mag- 
nitude of  the  error  involved.  An  erosive  burning  environment  was  simulated. 

It  was  necessary  to  have  a thin  sample  and  short  bum  time  due  to  the  finite 
volume  of  the  acoustic  tube.  The  sample  volume  was  the  same  as  the  strand 
volume  used  in  determining  the  direct  combustion  noise. 

The  sample  holder  was  constructed  to  hold  two  V'  x k"  x IV  strands 
(6,35  X 6,35  X 38,1mm)  with  a separation  of  up  to  V'  between  the  strands. 

This  holder  is  shown  in  Fig,  III  B-1,  Two  types  of  tests  were  conducted 
using  this  U-shaped  stainless  steel  sample  holder.  Atmospheric  tests  with 
plexiglas  sides  on  the  sample  holder  were  conducted  to  obtain  motion  pic- 
tures of  the  burning  samples,  Tliese  movies  allowed  surface  area  - time 
histories  to  be  obtained.  Stainless  steel  sides  were  used  when  the  samples 
were  burned  at  300  psi  in  the  acoustic  tube.  The  assembled  sample  holder 
for  the  pressure  tests  is  shown  in  Fig,  III  B-2. 

The  initial  assumption  was  to  design  the  U-shaped  channel  so  that  the 
exit  Mach  number  of  the  hot  gases  would  be  0.5,  if  the  sample  burned  normal 
to  the  V'  direction  uniformly.  It  was  hoped  that  by  controlling  the  amount 
of  ignitor  paste  the  initial  burning  surface  area  could  be  controlled.  Using 
available  thermochemical  equilibrium  data  for  the  UTP  3001  propellant,  the 
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Figure  III  B-1,  Apparatus  to  Test  Combustion-Flow  I 

Interaction. 
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initial  flow  area  was  calculated.  The  stagnation  temperature  of  5382°R 
gave  a speed  of  sound  of  4102  ft/sec  in  the  propellant  product  gases.  A 
realistic  burn  rate  from  strand  tests  was  ,35  in/sec.  In  order  to  obtain 
a Mach  number  of  0,5  at  the  exit  plane  of  the  U-shaped  channel,  the  initial 
sample  separation  of  the  two  strands  should  be  0,057"  (1.45mm),  This  initial 
channel  was  machined  to  allow  for  a V'  (6,35  mm)  separation.  This  large  sep- 
aration was  taken  up  with  a V'  * V x 1^'  (6.35  x 6,35  x 38,lmro)  steel  spacer, 
A series  of  six  atmospheric  tests  were  conducted  using  the  cross  flow 
hardware  to  obtain  motion  picutres  of  the  burning  profiles  through  the  plexi- 
glas  windows,  A non-aluminized  AP-HTPB  propellant,  T-48,  was  used  in  these 
tests.  The  presence  of  metal  would  obscure  the  motion  pictures.  The  first 
two  tests  were  unsucessful  because  of  extensive  time  delay  in  the  sample 
ignition  and  high  camera  speed.  In  the  third  test  the  amount  of  ignition 
paste  was  increased  and  the  camera  slowed  down  to  500  fps  from  800  fps. 

This  sample  started  burning  on  the  end  and  burned  for  45  seconds.  The  sample 
stopped  burning  V'  from  the  end  of  the  U-shaped  channel.  The  ignition  paste 
layer  appeared  to  remain  intact.  For  the  fourth  test  the  sample  separation 
was  eliminated  and  a minimum  amount  of  Ignition  paste  was  used.  Ihis  again 
resulted  in  an  end  burning  sample.  The  fifth  and  sixth  samples  were  success- 
fully burned.  The  scrands  were  tapered  slightly  to  give  an  exit  sep-'ration 
of  ,0625"  (1.59iiBn). Sample  frames  from  the  movies  of  these  tests  are  shown  in 
Fig,  III  B-3,  The  side  profile  of  the  fifth  sample  is  shown  in  Fig,  III  B-4, 
Assuming  that  the  sample  burned  uniformly  across  the  width  of  the  U-shaped 
channel,  the  area  as  a function  of  time  was  calculated.  This  is  plotted  in 
Fig,  III  B-5  along  with  the  associated  correction  factor  to  be  added  to  the 
acoustic  signal  obtained  from  burning  strands.  According  to  theory,  if  there 
is  no  cross-flow  effect,  the  amount  of  noise  should  scale  linearly  with  the 
burning  area. 
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Figure  III  B-4,  Time  History  of  Burning  Surface  Positions 
for  T-48  at  Atmospheric  Pressure, 
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Figure  III  B-5,  Time  His 


Direct  combustion  noise  was  then  measured  at  300  psi  using  the  MC-170 
and  T-48  propellant.  The  U-shaped  sample  holder  with  stainless  steel  sides 
was  used  for  these  tests.  The  samples  were  prepared  in  exactly  the  same 
manner  as  for  the  two  successful  movies.  The  resultant  sound  pressure  levels 
are  listed  in  Table  1 " B-1  along  with  the  corresponding  strand  results.  The 
samples  in  the  channel  burned  slower  than  the  corresponding  strands.  This 
could  be  due  to  a quenching  or  heat  transfer  effect  of  the  stainless  steel 
walls  as  opposed  to  the  silicon  rubber  inhibitor  on  the  strand  samples,  llie 
samples  with  cross-flow  generated  more  direct  noise  than  the  strand  as  was 
expected,  due  in  part  to  the  larger  burning  surface.  Even  when  the  nominal 
correction  of  10.8  db  from  Fig.  Ill  B-5  is  considered  there  is  still  a 7 
to  7.5  db  increase  in  the  sound  output.  This  is  interesting,  but  it  is 
found  in  the  following  section  the  noise  is  still  too  low  to  be  of  im- 
portance in  the  vibration  problem. 

The  acoustic  signal  for  the  T-48  strand  and  U-shaped  channel  are  shown 
in  Figs.  Ill  B-6  and  7.  The  RMS  value  of  the  acoustic  signal  was  obtained 
by  playing  the  recorded  signal  back  through  a Hewelett  Packard  model  3400A 
RMS  voltmeter.  The  tape  recorder  speed  reduction  and  the  time  constant  of 
the  meter  gave  an  effective  time  constant  of  ,25  sec,  Ihe  sound  pressure 
level  was  consistently  decreasing  for  the  strand  results  in  Fig,  III  B-6, 

This  is  as  the  tube  pressure  and  temperature  were  Increasing,  For  the  sample 
with  cross-flow  the  sound  pressure  level  was  constant  over  most  of  the  test, 
or  increasing  slightly.  The  sound  pressure  level  for  both  tests  varies  within 
a band  of  2 to  3 db.  The  only  change  in  the  correction  due  to  area  fluctua- 
tions shown  in  Fig,  III  B-5  would  be  near  the  end  of  sample  burning  and  would 
still  be  to  small  to  make  a meaningful  shift  in  the  trend  of  the  sound  pres- 
sure history  to  agree  closer  with  the  strand  result.  It  must  also  be  kept  in 
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Figure  III  B-6,  Sound  Pressure  Level  History  for  Strand 
Test, 
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Figure  III  B-7,  Sound  Pressure  Level  History  ^or  U-shaped 
Channel  Test. 


mind  that  if  the  cross  flow  samples  burned  as  fast  as  the  strands  there 
would  be  an  even  larger  difference  in  the  sound  pressure  level  than  7 db. 

Since  the  burn  rate  was  down  by  « factor  of  two  over  what  was  assummed 
in  the  design  of  the  crossflow  hardware,  the  Mach  number  and  Reynolds  num- 
ber based  on  the  strand  separation  at  the  end  of  the  Flow  channel  were  re- 
calculated, These  calculations  are  included  in  Table  III  B-2,  The  exit  Mach 
number  was  less  than  .26  for  all  conditions  while  Reynolds  numbers  indicated 
a definitely  laminar  flow. 

These  tests  indicate  that  perhaps  a significant  weakness  of  the  theory 
could  be  in  a neglect  of  cross-flow  effects  in  noise  generation  near  the 
propellant  surface.  Further  quantitative  work  in  this  area  is  required, 

C,  Temperature  fluctuations 

It  was  desirable  to  gain  an  estimate  of  the  level  and  spectral  content 
of  the  temperature  fluctuations  existing  in  the  products  of  combustion  of 
solid  propellants.  Such  an  estimate  is  directly  required  to  estimate  terms 
IIj^  and  III  in  Eqs,  (4).  This  was  a difficult  task  because  the  temperatures 
are  extremely  high  and  any  thermocouple  small  enough  to  respond  to  any  small 
scale  fluctuation  would  be  susceptible  to  condensed  phase  product  impact. 
Previous  tests  had  indicated  that  temperature  measurements  could  be  made 
by  burning  the  solid  propellants  in  a plexiglas  tube  and  using  a thermo- 
couple located  near  the  wall,  where  it  would  measure  the  temperature  of  a 
cooled  mixture  of  product  gases  and  wall  ablation  product  gases  [Ref,  (11)]. 

The  initial  test  configuration  consisted  of  a ,5"  (1,27  cm)  inside  diam- 
eter plexiglas  tube  with  a ,125"  (3,2  ram)  wall  thickness.  Hie  initial  tests 
were  to  be  conducted  at  atmospheric  pressure  with  a non-aliuninized  propel- 
lant, It  was  decided  to  use  a Chromel-Alumel  thermocouple  that  was  commercially 
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available  in  a ,001"  diameter  wire.  This  type  of  thermocouple  had  a relatively 
fast  time  response,  l,e,,  a short  time  constant.  It  had  been  used  in  previous 
tests  to  measure  temperature  fluctuations  in  a can  combustor  [Ref,  (12)], 
Techniques  had  been  developed  to  correct  the  power  spectral  densities  obtained 
from  these  thermcouples  considering  the  time  constant  of  the  individual  thermo- 
couples [Ref,  (12)],  This  type  of  thermocouple  will  withstand  a temperature 
of  1372°C  (1605°k)  with  a thermoelectrical  output  of  54,875  millivolts  tdien 
the  reference  junction  is  maintained  at  0°C,  This  means  that  the  distance 
above  the  burning  surface  and  the  distance  from  the  wall  must  be  controlled 
to  prevent  thermocouple  destruction  in  every  test,  since  the  final  tempera- 
ture reached  in  adiabatic  combustion  would  be  on  the  order  of  2500*^  to 
3000°K.  A .097"  (2.46  mm)  ceramic  cylinder  was  placed  in  the  flowing  hot  gas 
to  deflect  any  condensed  material  away  from  the  thermocouple.  This  tempera- 
ture fluctuation  measuring  apparatus  is  shown  in  Fig,  III  C-l,  The  length  of 
the  ,5"  (1,27  cm)  diameter  propellant  sample  was  chosen  to  allow  the  thermo- 
couple to  reach  an  equilibrium  value,  A ,25"  (6,35  mm)  height  was  sufficient 
for  the  atmospheric  tests,  while  1"  (2,54  cm)  height  was  used  for  the  tests 
at  300,  600  and  900  psi,  A photograph  of  a nonaluminized  sample  is  shown  in 
Fig.  Ill  C-2, 

The  flow  pattern  around  a circular  cylinder  changes  as  a function  of 
ud-, 

Reynolds  number  (Re  = E.  ),  For  small  Reynolds  number  the  wake  is  laminar. 

As  the  Reynolds  number.  Increases  a regular  vortex  pattern  known  as  klrmAn 
vortex  streets  are  formed  (60<Re  <5000  ),  At  still  higher  Reynolds  number 
the  wake  becomes  turbulent,  [Ref,  (13)],  The  frequency,  f,  at  \diich  the 
vortices  are  shed  is  related  to  the  flow  velocity  and  cylinder  diameter  by 
the  Strouhal  number,  St  = -^  , Ihe  initial  assumptions  predicted  a frequency 
of  2000  Hz,  This  was  beyond  the  expected  response  range  of  the  thermocouple 
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Figure  HI  C-1.  Temperature  Fluctuation  Test  Apparatus 


Figure  III  C-2.  Photogri?ph  of  Temperature  Fluctuation 
Test  Apparatus. 
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and  no  further  consideration  was  given  to  this  frequency.  The  post  test 
calculations  have  yielded  some  different  values  which  will  be  noted  as  the 
tests  are  discussed.  The  major  asstunption  with  this  apparatus  is  that  die 
wake  filling  and  emptying  time  is  short  compared  with  a typical  period  of 
the  fluctuations.  Thus,  the  frequency  content  of  the  fluctuations  should 
be  small  compared  with  the  Strouhal  shedding  frequency. 

An  initial  series  of  fourteen  tests  were  completed  using  a non-alumin- 
ized  propellant,  T-48,  at  atmospheric  pressure,  llie  time  histories  of  the 
thermocouple  response  were  recorded  on  magnetic  tape  for  reproduction  and 
Fourier  analysis.  The  thermocouple  and  its  location  was  fixed  for  all  tests 
in  this  series.  The  tube  was  shortened  after  each  test.  The  first  test  that 
was  completed  was  for  a thermocouple  located  initially  3,9"  (9,9  cm)  above 
the  burning  surface.  This  thermocouple  registered  a maximum  average  tempera- 
ture of  525*k  , The  last  test  of  the  series  was  for  the  thermocouple  initially 
located  ,27"  above  the  burning  surface.  The  temperature  history  for  this  run 
is  reproduced  in  Fig,  III  C-3,  The  fluctuating  component  of  the  thermocouple 
response  is  also  shown.  This  was  obtained  by  passing  the  thermocouple  signal 
through  a high  pass  filter  (Krohn-Hite  model  3323  R ) set  to  pass  frequencies 
above  5 Hz,  The  upper  limit  imposed  by  the  tape  recorder  was  30,000  Hz,  The 
temperature  history  for  one  of  these  tests  completed  using  an  aluminized  pro- 
pellant, UTP-3001,  at  atmospheric  pressure  is  shown  in  Fig,  III  C-4, 

The  RMS  value  of  this  fluctuation  was  obtained  by  two  methods.  It  was 
passed  directly  through  a Hewlett  Packard  model  3400  A RMS  voltmeter  or  ob- 
tained by  integrating  the  temperature  spectral  levels  shown  in  Fig,  III  C-5, 
This  latter  method  was  preferred  because  the  thermocouple  respoi.se  could  be 
corrected  for  its  individual  time  constant,  as  outlined  in  Ref,  (,12).  The 
same  problem  that  plagued  earlier  measurements  of  direct  combustion  noise 
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Figure  III  C-3.  Temperature  History  for  a Non- 
aluminized  Propellant. 


reoccured  when  the  same  Krohn-Hite  filter  was  used  to  high  pass  tlie  tempera- 
ture fluctuation  signal.  The  60  Hz  signal  and  its  harmonics  were  obvious  in 
the  low  level  temperature  spectra.  They  were  reduced  to  a lower  level,  as 
before,  by  subtracting  background  noise  spectra.  This  correction  was  not 
as  effective  since  the  thermocouple  output  changed  over  a much  larger  range 
for  an  individual  test  than  the  pressure  signal.  The  temperature  spectra 
were  corrected  for  both  time  constant  and  60  Hz  noise  before  the  trace  was 
integrated  to  give  a representative  RMS  value  of  temperature  fluctuation,  T', 

The  geometry  of  the  test  set  up  for  T-48  non-aluminized  propellant  is 
indicated  in  Fig,  III  C-6  along  with  a representative  average  value  of  the 
absolute  temperature  of  the  gas  passing  over  the  thermocouple.  The  extremes 
of  the  fluctuations  have  been  indicated.  The  non-dimensionalized  RMS  tempera- 
ture fluctuation  T',  divided  by  the  average  temperature,  T,  is  also  shown 
as  a percent.  For  the  nonaluminized,  atmospheric  tests  this  fluctuation 
varied  between  .8  and  1,4%,  The  burn  rate  for  these  tests  varied  from  ,016 
to  ,05  in/sec  (.4-1,2  mm/sec).  This  was  consistent  with  the  burn  rates  ob- 
tained from  the  U-shaped  channel  atmospheric  tests  of  the  T-48  propellant. 

The  results  for  the  absolute  temperature,  it’s  extremes  and  the  per  cent 
fluctuation  for  the  aluminized  propellant,  uTP-3001,  are  shown  in  Fig,  III  C-7, 
This  fluctuation  varied  from  4,4  to  6,9%,  The  burn  rate  for  these  tests  varied 
from  ,05  to  ,06  in/sec  (1,3 -1,6  mm/sec),  A photograph  of  the  last  test  of 
UTP-3001  at  an  initial  distance  of  6,2”  above  the  burning  surface  is  shown 
in  Fig,  III  C-8, 

The  burn  rates  for  these  tests  were  lower  than  that  assumed  in  the  initial 
calculations  of  Strouhal  ninnber  and  vortex  shedding  frequency.  Considering  a 
burn  rate  of  ,05  in/sec  (1,2  ttm/sec)  at  atmospheric  pressure,  the  velocity  of 
the  product  gas  for  T-48  (non-aluminized  propellant)  was  60  ft/sec  (18,3  cm/sec). 
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Figure  III  C-6.  Temperature  Measured  above  the  Burning 
Surface  of  a Non-aluminized  Propellant. 
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Figure  HI  C-7,  Temperature  Measured  above  the  Burning 
Surface  of  an  Aluminized  Propellant, 


The  resultant  Reynolds  number  was  1,6*  This  indicated  a laminar  flow  well 
below  the  lower  limit  for  the  formation  of  vortices  in  the  wake  of  the  ceramic 
cylinder.  For  the  aluminized  propellant  (UTP-3001),  the  flow  velocity  was  up 
to  180  ft/sec  (58,9  m/sec).  The  Reynolds  number  was  4,1  and  the  flow  shown 
in  Fig,  III~8  should  still  be  laminar  with  no  vortices  in  the  wake  of  the 
cylinder. 

It  was  concluded  from  these  tests  that  the  temperature  fluctuations  in 
the  product  gases  of  the  aluminized  propellant  were  considerably  higher  than  in 
the  non-aluminized  product  gases.  This  was  expected  because  of  the  burning 
aluminum  agglomerates  in  the  product  gases.  It  was  decided  to  conduct  further 
test  of  the  aluminized  propellant  at  300  psi. 

The  pressurized  temperature  fluctuation  measurements  were  conducted  in  a 
nitrogen  environment  in  the  acoustic  tube  shown  in  Fig,  III  A-1,  The  sample 
was  prepared  and  placed  in  a ,5"  (1,27  cm)  inside  diameter  plexiglas  tube 
as  for  the  atmospheric  tests.  It  was  necessary  to  increase  the  sample  height 
to  1"  (2,54  cm)  to  compensate  for  the  increase  in  burn  rate  as  the  pressure 
increased.  This  was  necessary  to  obtain  a sufficient  temperature-time  history 
for  Fourier  analysis.  The  plexiglas  tube  axis  and  the  direction  of  the  sample 
bum  rate  were  now  horizontal  instead  of  vertical,  due  to  the  acoustic  tube 
mounting  system. 

Two  initial  tests  were  conducted  at  a distance  of  9,9"  (25,1cm)  away  from 
the  burning  surface.  Both  thermocouples  apparently  reached  an  equilibrium 
temperature  (1142°K  and  1385^)  before  failing.  This  was  substantially  below 
the  upper  limit  of  the  Chrome  1 /A lumel  thermocouple  (1645^),  It  was  believed 
that  failure  was  due  to  particle  impingement,  A new  plexiglas  tube  was  used 
for  the  remaining  test.  The  length  was  increased  to  allow  an  initial  separation 
of  24"  (60,96  cm)  from  the  thermocouple  to  the  burning  surface.  Successful 
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tests  were  obtained  at  24,  18  and  12”  (61,  45,7  and  30,5  cm)  The  results 
for  these  three  tests  are  shown  in  Fig,  III  C-9,  The  burn  rate  of  these 
samples  varied  from  ,34  in/sec  (,86  cm/sec)  to  1,0  in/sec  (2,54  cm/sec). 

The  slow  burning  sample  was  for  the  test  with  24”  (61, cm)  of  separation 
between  the  thermocouple  and  burning  surface,  Hiis  sample  definitely  had 
problems  reaching  a steady  state  bum  rate.  The  slow  start  of  bum  did  not 
effect  the  resultant  temperature  measurements.  The  temperature  fluctuation 
varied  from  0,5  to  2,7%  This  was  considerably  below  the  values  obtained  in 
the  open  atmosphere  tests  shotm  in  Fig,  III  C-7, 

One  atmospheric  test  was  conducted  using  this  same  test  configuration, 

Ihe  sample  was  burned  in  the  acoustic  tube  that  was  flushed  with  nitrogen 
before  the  sample  was  ignited,  Ihe  absolute  temperature  was  consistent 
with  earlier  results,  but  the  flucttiatlon  was  less  than  indicated  in  Figure 
III  C-7,  The  temperature  fluctuation  due  to  after-burning  of  aluminum  in  a 
product  gas  mixing  with  nitrogen  may  be  less  than  that  produced  when  the 
product  gas  mixes  with  air. 

The  electrical  treatment  of  the  thermocouple  signal  was  improved  over 
the  course  of  this  experiment.  All  thermocouple  cables  and  temperature  com- 
pensation devices  were  electrically  shielded.  Ground  loops  and  60  Hz  noise 
pickup  were  minimized  by  device  and  cable  placement.  The  importance  of  the 
above  precautions  did  not  become  known  until  this  last  phase  of  the  investi- 
gation, This  was  consistent  with  the  elimination  of  the  high  pass  Krohn-Hite 
model  3323R  electronic  filter  from  the  experiment. 

Two  more  tests  were  conducted  at  19”  away  from  the  burning  surface  at  600 
and  900  psi,  Hie  absolute  temperature  and  the  temperature  fluctuation  have 
been  included  in  Fig,  111  C-9,  The  temperature  spectra  of  the  900  psi  test 
indicated  discrete  frequencies  at  700  and  820  Hz  as  shown  in  Fig,  111  C-lO, 
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This  suggested  the  possibility  of  a vortex  frequency.  Calculations  were 
made  knowing  the  test  sample  bum  rates  and  the  results  are  summarized  in 
Table  III  C-1,  The  resultant  Reynolds  numbers  indicated  that  there  could 
be  discrete  frequencies  in  four  of  the  six  tests.  Hiey  were  observed  in 
the  three  tests  where  the  maximum  frequency  considered  in  the  Fourier 
analysis  was  greater  than  the  frequency  indicated  by  the  Strouhal  number 
calculation.  The  occurance  of  more  than  one  frequency  is  not  suprising 
since  the  temperature  spectra  represents  only  10  samples  taken  as  the 
sample  bums  over  0,67  seconds.  A change  in  velocity  of  + 15%  would 
account  for  the  different  discrete  frequencies  indicated  in  this  table. 

The  major  point,  however,  is  that  the  peak  in  the  temperature  fluctuation 
spectra  occur  at  frequencies  well  below  the  vortex  shedding  frequency.  Con- 
sequently, the  spectral  shape  and  overall  rms  magnitude  are  believed  rep- 
resentative of  what  would  be  seen  if  the  protective  cylinder  were  absent. 

As  will  be  seen  in  the  next  section,  a typical  distance  between,  say, 

4 

400  p,m  diameter  burning  particles  Is  2,4  x 10  p,m.  For  the  data  appropriate 
to  Fig,  III  C-10  the  flow  velocity  was  less  than  50  ft/sec  (less  than  because 
the  gases  have  been  cooled  by  the  tube  and  are  less  dense  than  the  calculated 
adiabatic  flame  value).  Thus,  an  upper  limit  on  the  burning  particle  passing 
frequency  Is  690  Hz,  Tliis  is  the  same  ballpark  as  the  region  of  maximum 
spectral  density  on  Fig,  III  C-10,  Moreover,  it  was  observed  in  the  experi- 
ment that  as  the  pressure  was  lowered  and  the  gas  velocities  consequently 
raised  the  frequency  content  shifted  to  higher  values.  These  results  are 
consistent  with,  but  do  not  prove,  that  the  temperature  fluctuations  are 
caused  by  convection  of  high  temperature  combustion  fields  past  the  thermo- 
couple and  these  combustion  fields  are  attached  to  the  droplets.  It  is  known 
that  the  burning  metals  carry  their  own  diffusion  flame  with  them  in  rocket 
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Table  III  C-1  Reynolds  Number  and  Vortex  Shedding 

Frequency  for  Temperature  Fluctuation  Tests 


Pressure 

(psi) 

Burn 

Rate 

(in/sec) 

Distance 

(in) 

Velocity 

(ft/sec) 

Re 

f 

(Hz) 

frequency  indicated 
from  Fourier  A. 

Atmospheric 

.50 

12 

1040 

475 

26,900 

300 

1.0 

12 

33 

15 

300 

1.0 

18 

101 

47 

1500 

1720 

300 

.34 

24 

101 

47 

1500 

600 

.94 

19 

52 

47 

766 

705,820,890 

900 

1.50 

19 

48 

67 

843 

700,820 
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environments,  so  that  this  source  of  the  temperature  fluctuations  at  a point 
is  a plausible  one.  This  is  the  view  that  will  be  adopted  in  the  next 
section  and  is  one  which  critically  affects  the  outcome  of  the  calculations. 

Another  point  to  notice  is  that  the  temperature  fluctuation  levels  were 
not  particularly  high.  Part  of  this  may  have  been  caused  by  the  fact  that  the 
ceramic  shield  prevented  full  realization  of  the  actual  T',  On  the  other 
hand,  thermochemical ly  the  temperature  carried  in  the  diffusion  flame  sur- 
rounding a metal  droplet  should  be  significantly  higher  than  that  of  the 
surroundings  so  that  higher  T'  values  should  have  been  seen.  The  explana- 
tion here,  which  is  consistent  with  the  current  explanation  of  the  T' source, 
is  that  there  is  such  a large  distance  between  droplets  compared  with  a 
droplet  (and  flame)  dimension  that  the  rms  value  of  T'  is  suppressed  over  its 
peak  value. 


IV.  Order  of  Magnitude  Calculations 


As  mentioned  in  Section  II,  the  ultimate  output  desired  is  the  mean 
square  (or  rms)  pressure  fluctuation  and  its  spectral  distribution  on  the 
chamber  walls.  The  formulation  in  the  Appendix  is  in  terms  of  the  density 
fluctuation.  This  may  be  converted  to  a pressure  fluctuation  by  multipli- 
cation of  c„^  in  regions  where  a is  zero,  as  it  is  assumed  at  the  head 
and  nozzle  ends.  By  Eq.  ( 2 ) of  Section  II,  the  quantity  is  desired 

and  it  has  already  been  mentioned  that  an  initial  assumption  of  zero  cor- 
relation between  the  various  source  terms  of  Eqs.  (38)  of  the  Appendix 
will  be  made.  Consequently,  except  for  numerical  factors,  the  G's  of 
Eq,  ( 3 ) of  Section  II  are  made  up  of  times  the  square  of  the  indi- 
vidual numerical  terms  of  Eqs,  (38)  of  the  Appendix.  Consider  the  follow- 
ing manipulations  of  each  of  these: 
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with  ~ if  Tj^j^  is  only  correlated  with  its  space-separated 

^11 ^^i  ^ over  a distance  which  is  short  compared  with  a typical  dis- 
tance over  which  the  Green's  function  varies,  the  above  may  be  approximated 
by 
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Ihe  second  integral  here  has  the  magnitude  of  the  squared  value  of  1^^  (y^) 

ou 

times  a volume  over  which  the  is  correlated,  which  itself  may  be  a 
function  of  the  frequency.  The  first  integration  merely  weights  the 
second  by  the  square  of  the  Green's  function  and  adds  up  all  the  correlated 
volumes.  Consequently,  the  order  of  magnitude  estimate  is 
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Similarly,  any  term  involving  an  area  integration,  say  tenn  is 


plane  \ ®eo- 


e c c 


exit  plane 


To  obtain  the  mean  square  value  of  pressure,  the  Fourier  transform  product 
is  multiplied  by  the  2tt  times  the  reciprocal  of  the  ensemble  average  time 


and  integrated  over  frequency 
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That  is. 


„ 2tt  p p * 

do, 


For  order  of  magnitude  estimation  purposes  the  Green's  function  x^ill 
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be  taken  in  the  low  frequency  analytical  form  of  Eq.  (5  )•  Die  motor 
parameters  will  be  taken  as  those  of  the  baseline  case  of  the  Ap- 
pendix,  with  m-  = 0,5.  The  crux  of  the  problem  is  estimation  of  the 
spectral  character  of  the  source  terms  a and 


A,  Turbulence  Noise 

There  have  been  no  direct  measurements  of  the  turbulence  intensity 

( 8 ) 

or  spectra  distributions  within  a solid  rocket  motor.  A simulation, 
however,  has  shovm  intensity  distributions  which  are  similar  to  those 
of  pipe  flow, ^ Consequently,  the  spectral  distributions  of  Ref.  (7  ) 
will  be  used  for  estimation  purposes.  First,  look  at  term  I and  the 
makeup  of 

(T,^)  = p 1 (u  +u')(u  + u')1 

^ 11  cu  *^e  L e ' ' e ^ Juj 

where,  in  accordance  with  high thill  theory,  the  primary  contribution 

to  the  fluctuations  in  Tj , an  those  due  to  , >e  velocity  fluctuations 

(the  density  fluctuations  are  accounted  for  in  the  combustion  noise 

term  and  it  may  be  shown  that  the  density  fluctuations  arising  in  T^^ 

are  much  smaller  than  those  in  term  II,  ),  The  Fourier  transform  of  the 

D 

fluctuating  part  of  then  becomes 

(T,,  ) = p 1 2 u u + (u'2-  < u'2  >)  1 

ll'u)  ^e  L e (1)  'ou  J 

This  consists  of  a linear  part  plus  a part  due  to  nonlinear  fluctuations 
in  velocity  The  first  term  is  the  analogue  of  "shear-noise"  and  the 
second  is  "selt-noise",  using  the  terms  of  jet  noise  theory, 
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The  quantity  needed  is 


(T, 


11 


) (T, , ) = p ^ |4u  ^ u 11*  + 2u  j (u  - <u  ^^  > )*  u 

u)  11  U)  ^e  (.  e uj  m e L U)  ou 

+ u * (u^-  < u'^  > ) 1 

tt)  U)  J 


+ (u^-  < u'^  > ) (u^-  < u'2  > ) *\ 

(U  U)  'J 


In  jet  noise  theoxy  it  has  beer,  estimated  that  the  self  noise  terms  are 
of  the  same  magnitude  as  rhe  first  term  above. The  reasons  lie  in 
the  facts  that  in  jets  a)  the  turbulence  level  is  higher  than  in  pipe 
flow  and  b)  there  are  many  more  of  the  nonlinear  self  noise  terms  in 
jet  theory  than  in  the  above  so  that  the  sheer  weight  of  numbers  comes 
in.  Here,  however,  there  are  only  two  terms  involving  the  non-linear  part 
and  mean  turbulence  levels  are  only  of  the  order  of  57.,  Consequently, 
the  first  term  will  dominate. 


(T,,)  (T,.) 

11  U)  11  (JU 


4 p u2 
He  e 


u u 
U)  U) 


For  estimation  purposes,  Laufers*  pipe  flow  spectra  for  a particular 
radial  location  are  shown  in  Fig,  IV  A-1,  Shown  are  both  the  radial  and  axial 
velocity  fluctuation  spectra,  with  the  axial  velocity  the  one  of  interest 
here.  It  is  known  that,  for  fully  developed  flow,  the  spectra  are  invari- 
ant with  pipe  diameter  and  flow  velocity  if  they  are  plotted  as  F(St)  vs 
St  = 2k  a.  Ihe  scale  is  shown  in  Fig,  IV  A-1  such  that 

CO 

< u'  ^ >/u  ^ = J F(St)  dSt 
0 

The  next  problem  in  evaluation  of  Eq,  ( 6 ) lies  in  determination  of 
V^  , If  Taylor's  hypothesis  is  followed,  that  the  turbulence  is  convected 
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by  the  mean  flow  speed,  the  frequencies  seen  by  an  observer  would  cor- 
respond to  a flow  velocity  divided  by  an  eddy  size.  The  correlation 

.3 

volume  would  be  therefore  proportional  to  uj.  However,  Taylor's  hypothe- 
sis must  fail  at  sufficiently  low  frequency  because  the  eddies  would 
have  to  be  of  infinite  size.  At  the  low  frequency  end  of  the  range  it 
appears  reasonable  to  choose  the  integral  scale  of  the  turbulence  as 
the  appropriate  length  scale,  A function  which  has  the  above  behavior 


V 

cor 


e 


and  has  the  transition  to  a behavior  like  uj  when  the  frequency  is 
reached  corresponding  to  convection  of  the  integral  scale  eddies  by 
the  mean  flow  speed.  From  Laufer's  pipe  flow  data  « 0,4 

Finally,  since  there  is  little  or  no  turbulence  at  the  head  end 
and  highly  developed  turbulence  at  the  exhaust  end,  the  full  volume 
V does  not  contribute  to  the  noise  in  Eq,  (6)  Arbitrarily,  V is  di- 
vided by  two  in  the  magnitude  estimate.  Then  the  operation  of  Eq.  (2) 
is  performed.  The  result  is 


<_^  > I 

2 ^a 
Pe 


This  result  shows  the  volume  distributed  turbulence  to  be  an  infini- 
tesimal noise  source 


Tlie  term  is  treated  in  substantially  the  same  manner  as  above. 
The  nozzle  admittance  condition  is  a^  = - M^(y-1)/2v  , appropriate  for 
a short  nozzle.  Here  only  a surface  area  integration  is  required  and 
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2 2 


(T, ,)  is  estimated  as  above  and  S = /(1+  'T"  St)  . The  result  is 

' 11  (U  cor  ® 2a 


that 


> I,  = 2 . 10-5 

Pc 


which  is  also  a small  percentage  fluctuation  in  pressure. 

Term  1^,  presents  a more  difficult  magnitude  estimation  problem. 
A quantity  required  is 


(Ti,)  = 

11  UJ 


Byi 


••2  • • I 

(u-  + U-  Uj  + u,  u ^ + u/u  - uTuT 
'1  li  1 lili 
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It  is  assumed  that  the  turbulence  motion  is  primarily  vortical  so  that 
first  term,  being  a divergence  of  the  vortical  velocity  fluctuation, 
is  zero.  The  second  term  is  actually  a sum  of  three  terms  involving 
the  nonlinear  self  noise  part  of  T^^,  The  third  term  is  the  linear 
part  and  consists  of  contributions  from  only  the  axial  velocity  fluc- 
tuations. When  one  then  estimates 


S * 


there  arise  sixteen  non-zero  terms  of  which  one  consists  of  the  linear 
part  squared,  six  come  from  the  product  of  the  linear  and  nonlinear  parts, 
and  nine  come  from  the  products  of  the  nonlinear  parts.  At  a five  percent 
turbulence  fluctuation  level  one  would  need  about  twenty  of  the  linear- 
nonlinear  product  terms  to  be  comparable  in  magnitude  to  the  square  of  the 
linear  parts.  One  would  need  about  400  of  the  nonlinear-squared  terms  to 
be  comparable  to  the  square  of  the  linear  parts.  Consequently,  only  the 
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square  of  the  linear  parts  is  retained  and 


2 2 
p'^ 

•^e  e 


dy  dy 


Now,  if  the  turbulence  is  convected  at  the  flow  speed  9u  /9y  ?»  iuju  /u  , 

U)  u)  e’ 

and  this  will  be  assumed  here.  The  magnitude  analysis  then  follows  the 
previous  approach  to  yield 


- 7 X 10- 


This  is  now  large  enough  to  be  a contender  for  noise-making  at  the  1%  level. 

It  has  been  independently  verified,  by  using  the  linearized  vorticity-acoustic 
approach to  this  problem,  that  term  is  the  major  term  arising  in  a vor- 
ticity-nozzle  interaction.  Thus,  two  different  approaches  yield  the  same  answer, 
giving  confidence  in  this  Lighthill  approach. 
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B,  Combustion  Noise 


Consider  now  the  term  whereby 


* 224  2/  0**  \ 

II  II  = k c ^ M ^ ( 57?  g dS  ) 

aa  *^6  e r\jw  u)  / 


w 


This  is  estimated  as  usual  as 


M ^ g g*  S 'S, 
r u)  tu  w 6 


where 


'S  - ~ 

^6  " S 


dS 
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dS  57?  57? 

U)  U) 


is  a quantity  related  to  the  quantity  measured  in  Section  III,  In  fact, 

'2',  = S G,  where  G was  the  source  quantity  measured,  A reasonable  fit  to  the 
6 P 6 t) 

UTP-3001  and  MC-170  data,  representing  the  maximum  noise  propellant,  is 

3.98  X lO"^  , . 

Q - sec /rad 

6 1 + ( jiL  y 
\ 42  J 


at  300  psia.  Again  using  the  low  frequency  form  of  the  Green's  function 
and  integrating  over  all  values  of  frequency  as 


< P > -7 

-.-y-  II  = 8 X 10 

P ^ ^ 


Thus,  combustion  noise  near  the  propellant  surface  is  a small  noise 
source.  Recall,  however,  that  crossflow  appears  to  increase  this  level. 

It  is  possible  that  this  could  become  significant  in  the  presence  of  cross- 
flow.  However,  considered  here  were  the  maximum  noise  propellants,  and,  as 
a consequence,  combustion  noise  generated  near  the  propellant  will  be  re- 
jected as  an  important  noise  source  in  general.  It  is  possible  that  this 
should  be  reexamined  in  the  future,  Ihe  scaling  rules  of  this  noise  should 
also  be  noted  in  Kef,  (10):  It  is  possible  that  under  some  conditions 
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of  pressure  and  burn  rate  this  noise  may  become  important. 

Now  consider  the  noise  generated  in  the  chamber  volume  by  the  burning 
suspended  metal.  The  source  term  in  IIj^  is  D a^/Dt  which  is  the  time  deriva- 
tion of  the  entropy  fluctuations  seen  by  an  observer  moving  with  the  mean 
fluid  speed.  By  the  second  law  of  thermodynamics  this  is  composed  of  two 
parts  - the  heat  release  rate  fluctuation  and  convection  of  mean  entropy 
gradients  by  the  observer  by  the  turbulence  velocity  fluctuations.  It  has 
been  assumed  in  the  treatment  that  mean  entropy  gradients  are  absent  so 
one  is  left  with  the  heat  release  rate  fluctuations.  One  will  see  a heat 
release  rate  fluctuation  moving  at  the  mean  fluid  speed  by  the  convection 
of  individual  metal  combustion  fields  past  the  observer  by  the  turbulence 
of  the  field  and  by  any  metal  droplet  motion  not  strictly  parallel  to  and 
with  the  same  velocity  as  the  mean  velocity.  It  is  presumed  that  each  metal 

droplet  burns  with  its  own  vapor  phase  diffusion  flame  wrapped  about  the 
(14 ) 

droplet,'  Consequently,  the  volume  associated  with  the  combustion/diffusion 
field  surrounding  the  droplet  is  of  the  order  of  the  droplet  volume. 

It  is  instructive  here  to  compare  magnitudes, The  number  of  droplets  per 
unit  volume  is  given  approximately  by 


"d  = 


P ^ 
^e  d 


TT  d 


V 


d Ps 


so  that  d^^^  is  a typical  linear  dimension  between  particles.  Using  numbers 
for  aluminum  metal  at  300  psi  and  assuming  an  aluminum  loading  of  167, 


d 


inf'^^d 


28 


Tlierefore,  droplets  are  well  separated  and  can  burn  individually.  The  ratio 
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to  a typical  eddy  size  is  given  by 


int  d 

1 1 


for  d^=100  p,m  and  a = 1 ft.  Therefore,  there  are  about  10  droplets  per 
eddy.  The  major  point  here  is  that  the  volumes  over  which  combustion  events 
are  correlated  are  extremely  small  compared  with  the  correlation  volumes  of 
the  macroscopic  turbulence. 

In  the  magnitude  estimate  (Da/Dt)  will  be  estimated  as  u)CT  . Then 

A * 9 •2  4 xp*  \ 

c ^ II,  II  Y'^  P 8 g I i cr  ) I ) 

e b b ' U)  (U  \ Ju)  /\Ju)  / 


A slightly  different  problem  exists  here  in  estimation  of  the  volume  integrals 
since  such  a large  distance  exists  between  correlated  volumes,  as  shown  in 
the  magnitude  estimates  above.  The  meaning  of  the  double  volume  integration 
is  the  a c*  magnitude  times  the  correlation  volume  times  the  total  volume 

U)  U)  

over  which  correlated  volumes  exist.  Here 


' * * 3 ! d 

a dV  a dV  « a a d,  V ‘ 
UjJu)  UtfUjd  '“j,- 


where  the  last  factor  accounts  for  the  fact  that  the  correlated  volumes  are 
widely  disbursed.  Then  using  the  low  frequency  form  of  the  Green's  function 


and  the  usual  substitutions 


d/  1 (2n)^  2 Wd'Olnf)' 

s u)  a a * 

Sg  (1.9-Yn)2 


Now,  unless  a spectral  shape  is  introduced  for  an  integration  is  im- 

possible to  perform  to  obtain  < p^^  > As  a first  approximation,  assume 


that  all  of  the  noise  is  concentrated  in  a narrow  frequency  band  near  ou^  and 
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assume  that  u)^  corresponds  to  roughly  400  Hz,  where  the  temperature  fluctua- 
tions  peak  for  a run  in  Section  HI*  Then  even  if  < a > is  as  large  as 
unity,  the  result  is 


II.  = 2 


X 10  ^ d , 


with  the  units  of  d^  in  feet.  It  is  clearly  seen  that  for  any  realistic 
droplet  size  < p*  is  negligible.  This  is  a direct  consequence  of  the 

fact  that  the  correlation  volumes  are  so  small;  there  is  an  extreme  amount 


of  cancellation  of  sound.  None  of  the  assumptions  in  the  magnitude  analysis, 


if  relaxed,  would  change  this  conclusion. 


G,  Entropy  Noise 

Now  consider  term  III  with  b^  = M^/2, 
assumption.  Here 

* 4 72-2 

III  III  Cg  « / k p 


consistent  with  the  short  nozzle 


CT  ds  f a*  dS 
J CU  (0 


Again,  there  is  a problem  here  in  the  double  area  integral  because  the  cor- 
relation' areas  are  disbursed  because  the  droplets  are  widely  separated.  It 
is  here  being  assumed  that  the  temperature  fluctuations  behave  in  accordance 
with  the  experimental  results  of  Section  III  and  that  in  an  actual  rocket 
motor,  even  if  combustion  is  complete  by  the  time  the  exit  plane  is  reached, 
the  correlation  sizes  correspond  roughly  to  the  initial  metal  droplet  sizes. 
Then  , 
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Using  the  low  frequency  form  of  the  Green's  function  an  integration  may  be 
directly  performed  and  even  if  <^'2>j^sas  large  as  unity 


< 


HI 
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with  again  in  feet.  Clearly  entropy  noise  is  not  a factor  for  any 
realistic  d^. 

The  conclusion  is  evident  that  only  one  noise  source  is  clearly 
dominant  - that  due  to  vorticity-nozzle  interaction.  It  should  be  re- 
marked that  the  distributed  turbulence  is  close  in  magnitude,  but  upon 

/2  4 

close  examination  <.p  X “ M , Since  a high  M case  was  purpose- 

^a  ^b  ® ® 

fully  chosen  for  this  case,  a higher  number  than  normal  has  appeared  for 
the  distributed  turbulence  estimate.  In  addition  there  is  some  concern 
about  the  estimate  for  combustion  noise  near  the  propellant  surface  - if 
crossflow  alters  the  results  significantly.  Consequently,  comparison  of 
the  results  with  a motor  firing  are  desi’*eable. 
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V,  Motor  Prediction 


The  calculation  of  predicted  pressure  fluctuations  in  the  solid  rocket 
motor  interior  will  now  proceed  more  accurately  on  the  basis  of  only  one 
noise  source  - that  of  vorticity-nozzle  interaction.  To  obtain  a more  accu- 
rate calculation  than  that  produced  in  the  order  of  magnitude  calculations 
several  things  will  be  done.  These  are  a)  a more  accurate  spectral  distri- 
bution and  intensity  of  turbulence  will  be  produced,  b)  the  Green's  function 
will  be  accurately  calculated  and  c)  the  effects  of  propellant  feedback  re- 
sponse will  be  included. 

The  primary  interest  is  at  the  head  end  of  the  motor.  The  "exact"  formula 
for  the  spectral  density  of  the  pressure  there,  assuming  only  vorticity-nozzle 
interaction  noise  is  present,  is 


r - 2 2 2 
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In  order  to  evaluate  this  exactly,  the  temporal  and  spatial  behavior  of  the 
Uj^  component  of  turbulence  is  required.  In  Laufer's  pipe  flow  data  the  spec- 
tra of  Uj^  are  given  as  a function  of  space  location,  but  the  spatial  cross 
correlations  are  not.  In  the  study  of  Ref,  (8)  only  intensity  distributions 
are  given,  with  no  consideration  of  spectra  or  spatial  cross  correlations. 

It  is  consequently  necessary  to  introduce  into  Eq.  (7)  one  of  the  approxima- 
tions used  in  Section  IV  - the  one  dealing  with  an  estimate  of  correlation 
area  behavior  with  frequency  and  the  assumption  that  the  Green's  function 
is  invariant  over  a correlation  volume.  Equation  (7)  becomes 

* 

(8) 


G = — u:^  g g*  S S u u' 
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where  the  correlation  volume  will  be  estimated  as  before  with  S 


cor 
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The  spectral  behavior  of  will  be  taken  from  Laufer's  pipe  flow  data 

but  corrected  in  level  to  yield  the  intensity  distribution  of  Ref,  (8),  First, 
Fig,  V-1  displays  the  pipe  flow  spectra  data  as  a function  of  radial  position 
with  wavenumber  as  a parameter.  The  cross-section  is  broken  up  into  five  equal 
area  sections  with  midpoint  r/a  values  equal  to  0,22,  0,54,  0,71,  0,84  and 
0,95,  Values  of  the  spectra  from  Fig,  V-1  at  these  points  and  for  the  values 
of  St  on  Fig,  V are  stored  for  computation.  The  units  of  F are  immaterial 
since  only  shape  is  required  here.  Secondly,  the  spectral  absolute  level  at 
each  radial  location  is  scaled  such  that  the  area  under  the  curve  will  give 
the  proper  turbulence  intensity.  This  is  done  by  observing  that  the  experi- 
ments of  Ref,  (8)  yield  a maximum  intensity  near  the  wall 


, /2^.1/2  n 0/0  “0*8  / 

( < 'L,  >)  = 0,248  u m/s 

1 ■'max 


(9) 


with  u in  m/s.  It  is  also  found  that  a reasonable  correlation  of  the  mini- 
mum intensity,  at  r = 0,  is  given  by 


( < UL  = 0,37  (<  in  1 

T.  'min  L i max  J 
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again  with  ail  units  in  m/s.  It  is  next  observed  that  the  maximum  intensity 
occurs  so  close  to  the  wall  that  for  the  calculations  here  a monotonic  vari- 
ation from  the  minimum  to  the  maximum  will  be  assumed.  The  variation  which 
fits  the  daia  of  Ref,  (7)  is  a law  which  varies  the  intensity  according  to 
(r/a)^*^.  Thus  the  formula  for  radial  variation  of  intensity  is 
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Thus  at  each  radial  location  G^/lias  the  shape  given  by  Fig,  V-1  and  is 
scaled  so  that  the  area  under  the  spectral  plot  yields  Eq,  (11), 

The  next  problem  is  the  incorporation  of  a feedback  response  law  into 
the  Green’s  function  in  Eq,  (8),  A perusal  of  the  combustion  instability 
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lx  .erature^^^^ for  the  pressure  coupled  response  functions  (velocity 
coupling  IS  specifically  excluded  by  assumption)  yields  some  interesting 
observations.  First,  there  is  a paucity  of  data  in  the  low  frequency  re- 
gime of  interest  here  (say, <1000  Hz),  Secondly,  there  is  little  data  to 
support  differences  in  response  functions  between  metalized  and  non-metal- 
ized  propellants.  Thirdly,  there  is  no  adequate  theory  that  fits  the  data 

on  even  urumctali-jed  propellants.  Fourth,  a reasonably  significant  collapse 

2 

of  some  data  occurs  if  one  plots  the  value  of  p,  vs,  f/r  rather  than  f alone. 
This  last  f-'ct  suggests  that  the  thermal  wa\2  in  the  solid  phase  controls  at 
least  the  low  frequency  behavior  of  since  theories  based  on  such  a mechan- 
ism yield  an  f/r^  scaling  rule,^^^^  It  is  also  found  a posteriori  that  the 
feedback  ^ iiOoen  does  not  influence  the  results  greatly,  as  long  as  one 
is  sufficiently  removed  from  a stability  limit.  Consequently,  the  "A-B" 
model  of  Ref,  (16)  is  chosen.  This  sets 

t'  _ AB 

n * \ + A/X  - (1  + A)  + AB 
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with  A and  B free  parameters.  Best  fits  to  much  data  yield  B of  order  unity 
and  A of  the  order  of  20,  For  the  computations  to  follow  the  values  B = 1, 

A = 15  have  been  used  which  are  appropriate  for  A-35  propellant,''  ^ It  is 
important  to  notice  that  Eqs,  (12)  yield  the  physically  expected  result  that 
^ n as  u)  - 0,  The  calculation  of  the  Green’s  function  follows  from  the  formu- 
las in  the  Appendix. 

Before  proceeding  with  the  actual  calculations  it  is  first  important  to 
notice  a general  result  from  Eq,  (8),  Since  1 /a  Is  expected  to  be  relatively 
invariant  with  motor  design  varicbles,  the  maximum  relative  intensity  of  tur- 
bulence is  relatively  iiv^ensitive  to  u in  Eq,  (9),  and  noticing  the  approximate 
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low  frequency  form  of  the  Green's  function  in  Eq,  (5),  it  follows  that  an 
integration  over  frequency  of  Eq,  (8)  to  yield  <p^  >will  ”-eld  <p'  >/p 
completely  insensitive  to  motor  design  variables.  Therefore,  the  calcula- 
tions made  here  for  a particular  motor  will  be  more  or  less  general  for 
the  relative  pressure  fluctuations.  The  spectral  distribution  will  depend 
upon  motor  variables  but  < p >/p  will  not. 

It  should  be  noted  that  the  numbers  to  be  calculated  will,  of  course, 
depend  upon  the  input  intensity  of  turbulence.  Considered  here  are  motors 
with  a relatively  constant  cross  section  area  with  axial  position.  If  a 
complex  grain  design  yields  turbulence  levels  at  variance  with  the  inputs 
then  different  pressure  fluctuation  levels  will  be  achieved.  Using  the  phys- 
ical notion  that  complex  geometries  will  increase  the  turbulence  level,  it 
appears  plausible  that  what  is  being  calculated  is  a lower  limit  to  the 
pressure  fluctuations. 

Motor  data  were  supplied  by  Hercules  Incorporated/Allegany  Ballistics 
Laboratory,  A cutaway  view  of  the  motor  is  shovm  in  Fig,  V-2,  The  motor  con- 
forms nicely  to  the  assumptions  of  the  theory  in  that  the  great  majority  of  the 
grain  length  is  a center  p rtorated  grain.  The  motor  and  firing  designation 
was  LCLM-7,  The  propellant  formulation  and  characteristics  are  as  follows 


CTPB 

13,0% 

Fe203 

1,0% 

A1  (6  ^m) 

6,0% 

AP  (400  31 X,  200  p,m 

19X,  50  u,m  30%) 

80,0% 

* 

c 

5109  ft/sec 

r @ 1000  psi  70°F 

0,594  in/sec 

n 

0,406 

The  data  of  interest  were  over  the  first  0,8  sec  of  firing  vdierein  the 
motor  pressure  is  roughly  constant  at  1900  psia.  Other  numbers  required 
for  computation  are  as  follows: 
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c = 3630  ft/sec  M = 0.1447 

e e 

1 = 2.96  ft  a = 0.097  ft 

M = 0.00244 
r 

The  thermal  diffusivity  of  the  propellant  was  taken  as  the  value  for  pure 
-5  2 

AP,  Q-p  = 6 X 10  in  /sec.  Computed  values  for  the  maximum  and  minimum 
relative  turbulence  intensity  were  9%  and  0.87.,  respectively.  The  short 
nozzle  approximation  was  used,  which  is  perfectly  valid  for  the  frequency 
range  covered.  As  mentioned,  the  feedback  parameters  were  taken  as  B = 1,0 
and  A = 15,  This  is  a low  metal  loading  motor,  and  in  view  of  the  results 
in  the  appendix  that  the  metal  damping  affects  the  results  very  little, 

Yj  = = 0 are  chosen. 

The  results  are  shown  in  Fig,V-3  for  the  theoretical  spectrum.  The  results 
are  presented  as  spectral  level/50  Hz  bandwidth  to  be  consistent  with  the 
later  experimental  results.  The  spectrum  has  been  normalized  by  the  square 
of  the  mean  pressure  so  that  the  area  under  the  curve  yields  the  square  of 
the  fractional  pressure  fluctuation,  Tlie  number  for  the  rms  pressure  fluctua- 
tion is  p^  rms/P  ~ 0,0036,  about  the  same  level  as  estimated  in  Section  IV. 
Most  of  the  pressure  fluctuation  comes  from  the  spectral  distribution  at 
very  low  frequency.  The  first  longitudinal  mode  hump  contains  only  a small 
amount  of  the  mean  square  pressure, 

Tlie  experimental  spectrum  was  jwupplied  on  Ll»e  basis  of  200  points  spaced 
2.5  Hz  apart,  from  0-500  Hz,  The  pressure  pickup  was  a water  cooled  Kistler 
gauge  mounted  in  the  head  end,  lliere  was  no  detectable  signal  above  350  Hz, 
The  setup  of  the  electronics  was  actually  to  monitor  a + 100  psi  oscillation 
so  that  all  signal  here  is  very  close  to  the  noise  level  of  the  equipment. 
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FREQUENCY,  HZ. 

Figure  V-3.  Comparison  of  Theoretical  and  Experimental 
Spectra  for  Hercules  Motor. 


However,  the  comparison  with  theory  is  highly  favorable  so  that  it  is  be- 
lieved the  experimental  spectrum  is  truly  indicative  of  the  gasdynamic  noise. 

To  smooth  the  spectrum,  twenty  adjacent  bands  are  averaged  so  that  in  Fig,  V-3 
an  experimental  spectrum  with  a 50  Hz  bandwidth  is  presented.  Given  all  of  the 
massaging  of  the  data  and  the  observations  above,  the  error  in  the  experimental 
spectrum  is  probably  of  the  order  of  a factor  of  two,  or  + 3dB,  Given  all  of 
the  approximations  of  the  theory  the  same  kind  of  error  may  be  expected  in 
the  theoretical  curve.  In  this  light,  the  agreement  between  theory  and  ex- 
periment on  Fig,  V-3  is  considered  remarkable. 

It  is  noted  that  the  rms  pressure  levels  agree  within  a factor  of  1,67, 

The  spectral  shape  at  low  frequency  is  well  predicted.  The  major  discrepancy 
is  beyond  300  Hz  where  the  experiment  does  not  show  the  cut-in  of  the  first 
longitudinal  mode.  The  significance  of  this  discrepancy  is  believed  to  be 
that  the  turbulence  spectrum  of  pipe  flow  does  not  strictly  apply  to  the 
rocket  motor  system.  In  fact,  with  mass  addition  from  the  propellant  surface 
one  would  expect  a thickening  of  the  mixing  layer  between  the  core  flow  and 
the  injected  flow  as  compared  with  the  transition  layer  between  core  flow  and 
zero  wall  flow  in  the  case  of  pipe  flow.  This  increase  in  mixing  length  scale 
should  be  reflected  in  greater  low  frequency  content,  as  compared  with  pipe 
flow  results.  An  interesting  point  here  is  that  it  may  be  possible  to  infer 
Che  turbulence  spectrum  from  the  pressure  fluctuation  results. 

The  fact  that  the  theory  somewhat  underpredicts  the  level  of  the  overall 
fluctuation  should  probably  not  be  given  too  much  attention.  It  would  be 
tempting  to  blame  it  upon  the  neglect  of  the  other  noise  sources.  However, 
as  mentioned  above,  the  theory  and  the  experiment  contain  inaccuracies  which 
are  probably  larger  than  the  observed  difference  in  pressure.  It  should  be 
recalled,  however  that  this  is  considered  a lower  limit  analysis. 
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VI«  Conclusions 


1,  There  is  a single  noise  source  which  is  dominant  in  producting  pressure 
fluctuations  within  a solid  rocket  motor  cavity.  This  noise  is  caused 
by  turbulence  encountering  the  exhaust  nozzle  and  is  dominant  vdiether 
or  not  the  propellant  is  metalized  or  is  composite.  The  only  caveats 
are  that  a)  there  is  a possibility  that  at  low  port-to-throat  the 
turbulence  distributed  through  the  chamber  volume  may  become  important 
as  a noise  source  and  b)  there  is  still  an  unknown  effect  of  turbulence 
interaction  with  the  flame  zone  near  the  propellant  surface,  which 

may  be  an  additional  noise  source, 

2,  The  noise  floor  due  to  the  above  mechanism  is  at  a level  of  roughly 
0,47,  fluctuation  of  chamber  pressure,  is  low  frequency  in  nature  and 
can  be  larger  than  the  stated  level  if  the  motor  grain  design  induces 
a high  turbulence  level  in  the  chamber  gases, 

3,  The  noise  level  is  predictable  but  the  spectral  distribution  needs 
some  further  experimental  work,  because  the  theory  predicts  a slower 
high  frequency  roll-off  than  is  indicated  experimentally, 

4,  For  simple  motor  geometry  the  noise  considered  is  relatively  invariant 
with  motor  design  variables,  when  quoted  as  the  magnitude  of  the  frac- 
tional pressure  fluctuation, 

5,  A simple  theory  of  aeroacoustics,  based  upon  the  Lighthill  approach, 
has  recovered  all  of  the  expected  physics  of  the  noise  problem.  It 
allows  a framework  to  calculate  noise  magnitudes  excluded  here  if 
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future  experimental  work  shows  the  neglected  sources  should  be 
included# 

6#  As  opposed  to  results  of  instability  analysis,  the  effects  of 
metal  and  metal  oxide  damping  and  of  the  propellant  feedback  re- 
sponse functions  are  weak  in  noise  analysis,  as  long  as  one  is 
reasonably  well  removed  from  a stability  limit, 

7.  If  the  source  frequencies  are  primarily  belov?  that  of  the  first 
transverse  acoustic  mode  of  the  chamber  gases,  the  pressure  oscil- 
lations and,  hence,  the  motor  vibration  is  purely  axial,  except  for 
manufacturing  misalignments  and  small  local,  rather  than  global, 
pressure  fluctuations. 
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Appendix  - General  Theory 


Formulation-Differential  Equation 

There  are  several  possible  theoretical  formulations  for  the  prob- 
lem i»f  the  interior  aeroacoustics  of  a rocket  motor,  Tliey  are  all 
approtimate  methods  and  the  method  chosen  should  represent  a balance 
between  analytical  simplicity  and  the  accuracy  of  the  data  to  be  used 
in  the  theory.  After  a review  of  several  possible  approaches  it  was  de- 
cided to  attack  the  problem  from  the  Lighthill  theory.  This  was  done  be- 
cause a)  an  analytical  solution  to  the  problem  was  anticipated  by  this 
method  b)  the  theory  is  known  to  give  proper  orders  of  magnitude  and 
scaling  rules  in  calculations  of  jet  noise^^^  and  combustion  noise, 
c)  the  theory  is  especially  good  for  relatively  low  Mach  number  flows 
such  as  found  in  the  interior  of  rocket  motors  d)  the  data  base  to  be 
used  in  the  theory  will  contain  inaccuracies  that  make  use  of  a more 
accurate  theory  unjustified  and  e)  it  was  found  possible  to  include  all 
of  the  anticipated  physical  effects  of  the  problem  through  this  approach. 

Writing  the  continuity  and  momentum  equations  for  the  gas  phase 


dt  ^ 9x.  ^P^i^ 


m 


(1) 


9t  j 9x^  i 


- ^ +f 
Sk.  i 


(2) 


85 


In  Eq,  (1)  the  mass  source  for  the  gas  phase  comes  from  the  combustion 

of  suspended  condensed  phase  materialo  In  Eq,  (2)  the  force  term  comes 

from  drag  of  the  condensed  phase  material.  Equation  (2)  considers  an  in- 

viscid  gas,  except  for  the  drag  term,  since  viscosity  plays  little  role 

(4) 

in  noise  generation;  the  larger  turbulent  eddies  are  responsible  for 
noise,  whereas  viscosity  is  primarily  active  in  the  fine  scale  structure 
of  the  turbulence, 

Ihe  Lighthill  equation  is  obtained  by  a combination  of  Eqs,  (1)  and 


(2).  This  is 


o -.2  ^2  ^2.  2 , 

2 e ax  ax  ax. ax  i y ax. ax, 

0^  X X XX  XX 


Equation  (3)  holds  for  any  constant  c^,  but  will  be  now  identified 
with  the  speed  of  sound  at  the  nozzle  entrance  plane.  The  last  two 
terms  involving  mass  addition  are  typically  dominated  by  heat  addi- 
tion, which  will  arise  from  the  third  term  on  the  right  hand  side  of 
Eq,  (3);  consequently,  the  last  two  terms  of  Eq,  (3)  are  dropped.  Since 
it  will  be  found  that  it  is  most  convenient  to  work  in  the  frequency 


domain,  rather  than  the  time  domain,  the  Fourier  transform  is  taken  of  Eq,  (3), 
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The  result  is 


9 p o 

r — ^ + k p 


(T.  ^ . _ PUi 

Bx.Bx  ^Ju)  Bx.Bx.  ^u)  ” 2 

1 j 1 i Cg 


+ 


Bx. 


(fl,> 


U) 


(4) 


In  Eq,  (4)  all  quantities  are  now  the  transform  of  the  fluctuating  quanti- 
ties; the  d.c.  components  have  been  removed  by  subtraction. 

Now,  two  kinds  of  condensed  phase  particles  will  be  considered  in 
the  calculation  of  f . , a metal  and  its  oxide.  It  is  furthermore  assumed 
that  a single  representative  size  may  be  assigned  to  each  particle  type. 
Since  the  metal  will  disappear  as  the  nozzle  is  approached,  due  to  com- 
bustion, the  size  assigned  to  the  metal  will  be  some  average  size  dur- 
ing the  metal's  lifetime.  An  average  number  density  will  also  be  assigned 
to  the  metal  and  oxide,  recognizing  that  in  reality  this  will  be  a function 
of  chamber  postion. 

If  the  number  density  were  a constant,  the  particles  were  suspended 
in  a stream  flowing  at  constant  velocity  u^,  and  the  particles  were  also 
flowing  at  velocity  u^,  there  would  be  no  drag.  Considering  fluctuations 
about  such  a condition 

ff  = n D.' 

1 1 

(5) 
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The  condensed  phase  momentum  equation  is 


4:  '’pi  ■•"i 


which  becomes  for  fixed  mass  m and  the  approximation  of  Eq,  (5) 


r— — + u ^ Vp  = - of 
dt  e dx  1 


Using  Stokes  law  for  the  drag  coefficient 


C 

D Re 


and  taking  the  Fourier  transform  of  Eq,  (6), 


iiu  + iai  » (V,  - V ) 

Pi^  Pi,., 


2 a (v.  - V ) 

lU)  p,- 


The  exact  solution  to  Eq.  (7),  subject  to  the  condition  that  there  is 
no  particle  motion  fluctuation  at  the  head  end  (x=0)  is 


-Bx  1*^  Bx  , 

V = a e V.  e*^  dx 

^im  o UJ 


with  B - iu)/up  + O'.  For  large  g,  an  asymptotic  solution  to  Eq.  (8)  is 


V = V a/B 
iu)  U) 


Equation  (9)  Is  exact  in  two  limits;  these  are  g "*  ® (either  a or  (u/u^ 
large)  or  a — o«  There  is  error  committed  for  moderate  g (both  a and 
lu/u^  moderate).  Because  of  the  exact  behavior  at  two  limits  this  error 
involved  in  use  of  Eq.  (9)  is  expected  to  be  minimal,  especially  in 
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view  of  the  other  approximations  leading  to  Eq,  (7),  Moreover,  it  will 
be  found  that  there  is  a relatively  small  effect  of  condensed  phase 
drag  on  the  solution.  Consequently,  Eq,  (9)  is  accepted  and  the  Fourier 
transform  of  Eq,  (5)  becomes 


f.  =nau(^-l)v. 

1 e ' B 1 

UJ  to 


(10) 


Neglecting  the  mass  addition  term  in  Eq.  (1),  consistent  with  a 
previous  assumption,  and  assuming  small  fluctuations  about  the  exit 
plane  conditions  together  with  a mean  one-dimensional  flow  of  the 
chamber  gases,  the  transform  of  Eq,  (1)  becomes 


1 to 


= 0 


whereby  with  the  aid  of  Eq,  (10) 


IJi'io.  ■ “ to*  ] i ) L ] 


■a  ^ ] 


(11) 


Finally,  a factor  of  (2iku^u)/c^)  Bp^/Sx  is  subtracted  from  both 
sides  of  Eq,  (4)  and  Eq,  (8)  is  substituted  in  Eq,  (4)  with  tae  addi- 
tional approximation  in  Eq,  (8)  that  u = u^.  The  result  is 


S p, 


2 


^x.  ■ "‘'^e  ^ [2+2d]  +K,^  Pu,  = " ^ 


r H 


—2  ~-r  (T^. ) + ( p - -^  ) - z i k M ^ 


= k^  +*d^ 


(12) 
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Some  of  the  approximations  in  the  first  of  Eqs,  (9)  can  now  be  ration- 
alized, The  operator  on  the  left  hand  side  is  easily  recognized  as  that 

appropriate  to  wave  propagation  in  the  presence  of  a uniform  mean  flow 

2 

where  a)  the  Mach  number  is  low  so  that  terms  of  the  order  of  M com- 

e 

pared  with  unity  have  been  neglected  and  b)  a correction  for  droplet 
drag  is  present.  The  insertion  of  for  the  actual  variable  M(x)  or 
M(x^)  is  an  approximation  which  is  valid  in  the  spirit  of  this  treat- 
ment, which  is  basically  an  order  of  magnitude  treatment.  The  main 
effect  of  the  second  term  on  the  left  hand  side  is  on  the  resonant 
frequencies  of  the  motor  (as  may  be  seen  later).  While  it  is  believed 
reasonable  to  incluae  this  effect,  it  is  done  in  onxy  an  approximate 

manner.  The  z,  term  in  the  second  term  on  the  left  hand  side  is  actu- 
d 

ally  small  compared  with  2.  As  may  be  deduced  by  direct  computation 
Zd  < which,  at  a maximum,  is  of  the  order  of  0,2,  Consequently, 
errors  in  the  treatment  of  condensed  phase  drag  affects  this  term  very 
little.  Droplet  drag  is  essential  in  the  last  term  on  the  left  hand 
side,  in  the  wave  number  K.  z^is a complex  number  and  the  imaginary 
part  will  induce  an  imaginary  part  to  K.  This  will  be  responsible  f r 
a damping  effect  of  condensed  phase  drag.  There  is,  in  fact,  an  opti- 
mum damping  condition.  This  is  seen  by  investigation  of  the  imaginary- 
part  of  z^  which  is  proportional  to 


au  /(i) 
e 


This  has  a maximum  at  a = Thus,  this  approximate  treatment  retains 

the  known  feature  of  more  exact  treatments  that  at  each  frequency  there 

(18) 

is  an  optimum  particle  size  for  damping.  Summarizing,  even  though 
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there  are  approximations  in  the  treatment  of  the  particle  damping  and 
flow  effects,  the  essential  physics  are  retained  and  the  scaling  proper- 
ties with  various  variables  are  retained. 


Equation  (9)  is  the  modified  Lighthill  equation  to  be  solved  for  this  in- 
terior aeroacoustics  problem.  In  the  spirit  of  that  approach,  the  right 
hand  side  is  considered  calculable  from  estimates  of  turbulence  behavior. 

The  right  hand  • ide  contains  the  Li^thill  quadrupoles,  Ty  the  combustion 
noise  term,^^  p^,  a pressure  fluctuation  term  which  will  be  addressed 

later  and  a term  involving  M . This  last  term  will  now  be  dropped  as  being 

e 

small  compared  with  the  other  source  terms  (right  hand  side  terms),  espe- 
2 

cially  the  9 P(ju term.  The  right  hand  side  of  Eq,  (12)  will  then 
act  as  the  forcing  function  for  the  wave  operator  on  the  left  hand  side, 
where  now  the  right  hand  side  is  set  equal  to  - F with 


r - 1.  (r  ^ ^ ^ ^ 

^ “ c 2 ax. ax.  ax.sx/P'"  c 2 ) 

I J 1 i e 


To  solve  the  equation  appropriate  boundary  conditions  are  required. 
Formulation  - Boundary  Conditions 

The  configuration  is  as  shown  in  Fig.  A-1 , The  differential  equation  in 
Eqs,  (12)  holds  inside  of  V,  which  is  enclosed  by  the  bounding  surface  S. 
At  this  point  it  will  be  assumed  that  the  cross-section  area  is  constant 
with  axial  distance,  , in  fact,  that  the  grain  is  cylindrically  per- 
forated, although  this  will  be  relaxed  later  to  allow  more  complex  cross 
sectional  shapes.  The  surface  at  the  head  end  coincides  with  the  wall, 
and,  since  flow  velocities  are  small  there,  it  is  assumed  there  is  no 
turbulence  there.  The  surface  bounds  the  nozzle  entrance  plane,  accross 
which,  of  course,  a turbulent  fluid  is  flowing.  The  surface  is  set  far 
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Figure  A-1,  Schematic  of  a Rocket  Motor  and  Nomenclature 
Illustration, 


enou^  off  of  the  propellant  surface  that  all  reactions  have  been  com- 
pleted except  those  of  any  metals.  Typically  this  distance  will  be  of 
the  order  of  hundreds  of  p,  m , Across  the  gas  flow  is  assumed  purely 
radial.  Moreover  it  is  assumed  that  the  turbulence  does  not  penetrate 
S^.  This  important  assumption  is  tantamount  to  assuming  that  there  is 
no  erosive  burning  effect, ^^^^The  assumption  is  made  for  simplicity  in 
the  analysis  to  follow  and  will  be  reexamined  when  some  cross  flow 
experiments  are  considered.  Across  S^,  therefore,  there  is  a radial, 
nearly  laminar  flow  containing  products  of  combustion  and  metals,  Ihe 
metal  will  be  considered  to  be  primarily  consumed  in  V,  Typically,  even 
for  heavily  metalized  propellants,  the  added  energy  by  combustion  of  the 
metal  is  less  than  20%  of  the  overall  energy  of  the  gases.  Consequently, 
in  V there  are  considered  gases  with  a mean  speed  of  sound  equal  to  c^ 
about  which  acoustic  and  turbulent  fluctuations  are  taking  place.  The 
gases  will  be  assumed  perfect  so  that 


Y 


e 


and  the  fluctuation  in  entropy  is  given  by 


c 

P 


(13) 


Equation  (12)  now  becomes 


dp  , 

- — ^ - (2  + ^)  ikM  + K p 

9Xi5x^  ^ e 9x 


(U 


’’e  dx  dx/^U) 
i i 


(14) 
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1 


where  the  Fourier  transform  of  Eq,  (13)  has  been  used. 

The  boundary  conditions  for  Eq,  (14)  must  be  developed  on 

and  S , On  S,  the  wall  is  impervious.  Assuming  f.  = 0 at  the  head  end 
e n 1 

because  velocities  and  drag  are  small,  Eq,  (2)  may  be  evaluated  in  the 
X-  direction  for  which  v^  = u = 0,  Consequently, 

tfkr.e-" 


At  the  nozzle  entrance  plane  a nozzle  admittance  condition  for  choked 
nozzles  may  be  used.^^^^It  will  be  found  later  that,  because  of  che  inter- 
est here  in  only  very  low  frequencies,  only  the  plane  wave  acoustic  mode 
will  be  under  consideration.  In  this  case  the  admittance  condi.toJis  is^^®^ 


b 

cu 


do 


= 0 


(16) 


The  oo  subscript  indicates  the  plane  wave  mode  and  is  equi. valent;  to  a 
cross  section  average  of  the  indicated  physical  quantity.  Considering 
a one-dimensional  mean  flow  at  the  nozzle  entrance  ami  co!iipleted  burning 
so  that  m = 0,  Eq.  (1)  may  be  written  as 


U)  p 


dx. 


U) 


+ u 


Bp 

Bx 


= 0 


0 ••) 


at  the  nozzle  entrance  plane  after  linearization  of  iq,  (1).  Because 
condensed  phase  drag  is  included  in  the  basic  differential  equation 
and  will  be  found  to  have  only  a minor  effect,  it  will  be  excluded  in 
writing  the  momentum  equation  at  the  nozzle  entrance  plane.  Linearizing 
Eq,  (2)  about  the  mean  flow  condition,  the  transform  of  Eq,  (2)  becomes 
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(18) 


9v. 


i «)  Pp  + Pp 
U) 


9p 

JiSi 

ax 


The  cross-section  average  of  Eqs.  (17)  and  (18)  can  be  taken,  yielding 

equations  in  u , p and  p , and  a combination  of  Eqs,  (16)-(18) 
^00  P*oo  ^0 

then  yields 


9p, 


oo 


ax 


e^UJoo 


q;  = -ik(M  +Ya) 
e e ' u) 


P = - p jr~°°  ^ ^ Pp  ^ (^  + ) 

^e  ^^e  dx  ^e  uj^^'  U)  ' uj 


(19) 


where  the  transform  of  Eq,  (13)  has  been  used  to  eliminate  p . Making 


U)i 


00 


the  Eq,  (13)  substitution  into  Eq,  (15),  taking  the  cross  section  average 
and  recognizing  that  = 0 at  the  head  end 


ap 

(U  \ 

9P)  . 0 

ax  Jo 


(20) 


The  boundary  condition  on  the  side  wall  requires  some  further  approxi- 
mations, Hie  fluctuation  in  gas  mass  flow  per  unit  area  crossing  may 
be  written  under  the  purely  radial  flow  approximation  as 


m 

Jli 

m 


(21) 


It  is  assumed  that  this  consists  of  the  sum  of  two  independent  contribu- 
tions - a fluctuation  due  to  the  propellant  heterogeneity  and  a part  due 
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to  a pressure  sensitivity  of  bum  rate.  Hiat  is 


m 

= 

m 


(22) 


The  feedback  response  factor,  |j,  , is  allowed  to  be  complex  and  frequency 

dependent.  However,  as  frequency  tends  toward  '^ero,  it  is  demanded  that 

p,  tend  toward  the  real  valued  pressure  exponent  of  the  burning  rate.  The 

factor  57^  will  be  determined  frcmi  experiment.  It  will  also  be  assumed 

that  a = 0 on  S so  that  the  oscillations  are  nearly  isentropic  on  this 
U)  w 

surface.  Neglecting  condensed  phase  effects  also,  the  transform  of  the 
fluctuation  of  the  radial  momenbrni  equation  may  be  written  as 


i k 


+ p 


M avj. 
r ^u) 

e Cg  hr 


1 5p  5p 

i , _Ui JJi 

c^^  ar  ar 


(23) 


The  second  term  of  Eq,  (23)  will  be  neglected  compared  with  the  first 
since  is  usually  small  (of  the  order  of  0,02)  Then  substituting  Eqs, 


(21)  and  (22)  into  Eq,  (23) 


9r  /x,a,G 


+ a 
w 


'y 


= i k pg 
= - i k p 


Mj.  (y’^T-l) 

M ^ 
e r u) 


(24) 


Although  it  may  seem  a bit  cavalier  to  throw  out  the  term  above  with 

the  multiplier  even  when  the  terms  in  Eqs.  (24)  are  retained  with 

the  same  multiplier,  the  reason  lies  in  the  essential  behavior  of  the 

term,  contains  the  feedback  term  and  cannot  be  discarded,  especially 

since  p,  can  be  large.  B contains  the  source  term  51^  for  v^ich  the 

w u) 
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magnitude  is  not  yet  known.  The  point  is  that  in  Eq.  (23)  the  second 
terra  is  expected  to  be  small  compared  with  the  first. 

Equations  (19)  ,(20)  and  (24)  are  the  boundary  conditions  on  the 
problem.  Although  Eqs.  (19)  and  (20)  are  specifically  written  for  the 
plane  wave  mode,  whereas  Eq.  (24)  is  general,  it  will  be  seen  that  this 
causes  no  difficulty.  For  solution  of  the  differential  equation,  therefore, 
the  boundary  condition  is  of  the  general  form 

(25) 


where  n^  = e^^ , e^  and  -e^  at  the  nozzle,  side  wall  and  head  end, respective- 
ly. 

Solution 

The  method  of  solution  chosen  is  through  use  of  a Green's  function. 

If  the  Green's  function  satisfies 


+(2«d)  IkM^ 


- y^) 


(26) 


with  the  boundary  conditions 


Sg 

_JU 

3x 

)o,r,0 

dg 

JlStt 

Sx 

)'t,r,Q 

dg 

—ill 

dr 

)x,a,0 

+ (2  + z^)  ikM^g^  = 0 


(27) 
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it  may  readily  be  shown  that  is  given  by 


p (x.)  = 1 g (x. ,y.)  r (y.)  dV(y.) 
(U  1-  J 0)  11  1 1 


+ g dS 
J tu 
S 


(28) 


lhat  is,  Eq.  (28)  is  the  solution  to  Eqs.  (12)  subject  to  Eq,  (25),  The 

solution  consists  of  a suposition  of  volume  distributed  sources  and  surface 

distributed  sources.  The  Green's  function  gives  the  effective  weight  of 

the  source  (eith.  r f or  at  the  observation  point  x^.  The  f tem  contains 

Tij  aed  a fluctuations  and  the  ^ terms  contain  a and  771  fluctuations.  Evi- 

denuly,  the  solution  for  g will  contain  the  nozzle  and  feedback  effects, 

(U 

as  well  as  those  of  condensed  phase  drag. 

The  solution  to  Eq,  (26)  subject  to  Eqs,  (27)  is  carried  out  by 

(21) 

standard  methods. 

The  solution  is 


00 

L. 


<1-)  (5.») 

m,n=0 
a=  + 1 


Y = j (K  5)  „,cp  a = + 1 
mn  m mn=  sin  m cp  a = "1 

mn  m inn  w m mn 

F " + (2  + z,)  ikM^  Fj  + (K^-K  ^ ) = 


mn 


e mn 


mn 


6(x-y)Y^°(r.O) 


r ' (0)  + (2  + zj)  ikM  F (0) 
mn  ' ' ' e mn 


(29) 
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v^ich  requires  the  solution  to  the  F differential  equation.  For  any 

tnn 

in,n  pair  the  solution  will  look  like 


The  lowest  root  of  the  defining  equation  for  the  in  Eqs,  (29)  is 

designated  as  the  00  mode  and  corresponds  to  the  nearly  plane  wave  mode. 

Noting  that  a is  small,  K w 0,  and,  by  expanding  the  Bessel  functions 
w oo 

in  their  appropriate  power  series  for  small  values  of  the  argument,  it 


is  found  that 


(31) 


which  is  a small,  complex  number.  The  next  highest  root  is  the  01  case 

which  corresponds  to  the  first  antisymmetric  transverse  mode.  Hie  solu- 

2 

tion  will  be  restricted  to  values  of  frequency  such  that  K < Kq^  • The 

transverse  modes  are  then  "cut-off"  because  k for  all  transverse  modes 

mn 

contains  an  imaginary  part  from  Eqs,  (30),  The  solution  for  F^^  then  shows 
that  these  modes  decay  with  x;  they  are  non-propagating.  Consequently, 
only  the  plane  wave  mode  need  be  considered,  which  is  the  first  term  in 
the  series  solution  for  g^,  in  Eqs,  (29),  In  the  following,  the  oo  sub- 
script will  be  dropped,  it  being  understood  that  this  is  the  only  mode 
under  consideration.  This  is  the  reason  for  expressing  the  boundary  con- 
ditions of  Eqs,  (19)  and  (20)  in  plane  wave  form. 
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Bie  solution  to  the  F equation  with  the  boundary  conditions  of  Eqs, 
(29)  also  requires  the  statement  that  F is  continuous  at  the  singular 


point  X = y.  Then 


F = A + A e^^ 

R “ + 

F = B e^-y  + B . e"^y 

Li  * T 


X > y 


X < y 


with  tiae  A's  and  B's  determined  from  the  four  simultaneous  equations 


B_(l-e)  - B^(H€)  --  0 


B (1+e)  - B^  e^^  - A_(l4e)  + A^(l-e)  0^"^= 


B e’'^-^  + B,  e\^  - A e"^-^  - A.  e'^’^  = 0 


A_(l+q)  . (l-q)  A eV  - 0 


e - (l+‘=j/2  ) k/k„„ 


Since  Y ss  1,  in  the  plane  wave  limit  g = F(x,y). 
oo  tu 

An  interesting  calculation  is  in  determining  g^  in  the  limit  of  low 
frequency.  Carefully  taking  the  limit  of  u)  “•  0 in  Eqs,  (33)  and  Eq,  (32), 
it  is  found  tint  a)  the  Green's  function  becomes  independent  of  x and  y 
(this  becomes  a bulk  mode  oscillation)  and  b)  the  following  simple  formula 


emerges : 


* 2 


“>  hJ  k^  ( ^ - yIT)^ 


2 2 

Thus,  since  g es  l/k  , the  Green's  function  becomes  unbounded  at  low 
U) 

frequency  which  physically  will  mean  that  low  frequencies  are  favored  by 
the  duct  acoustics.  Note  in  Eq,  (34),  if  y,  is  large  enough,  g^  -•  ®, 
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That  is,  with  enough  positive  feedback  from  the  propellant  the  system 
may  be  driven  unstable,  as  is  well  known.  Since  ~ may  actually  be  com- 
plex, the  square  operation  of  Eq,  (34)  is  really  to  be  interpreted  as 
the  quantity  in  brackets  times  its  complex  conjugate.  In  Eq,  (34)  the 
effect  of  droplet  drag  disappears  at  low  frequency.  This  low  frequency 
solution  will  be  found  to  be  useful  for  order  of  magnitude  arguments. 

In  general  however,  Eqs.  (33)  must  be  solved  exactly. 

Extension  to  Non-Circular  Geometry 

Since  it  has  been  found  that  only  the  plane  wave  mode  is  under  con- 
sideration, this  will  permit  an  approximate  generalization  of  the  results 
to  non-circular  shapes  of  (nearly)  constant  cross-section  area  as  a func- 
tion of  axial  position.  An  example  would  be  that  of  a star  gain.  If  an 
effective  circular  radius  is  defined  in  terms  of  the  cross-sec -ion  area 


it  is  assumed  that  a^  may  be  used  in  place  of  a in  the  theory  above.  The 
radius,  a,  appears  in  and  physically  measures  the  lateral  length  scale 
of  distortion  of  the  plane  wave.  The  same  physics  are  obtained  by  the  ap- 
proximate procedure. 

Assuming  for  approximate  calculation  purposes  that  there  are  no  con- 
densed phases,  a mass  balance  yields 


p u S = p V S 
•^e  e c '^e  r w 


(35) 


with  S = P-t,  This  yields  a constraint  on  the  ratio  of  u /v  j5s  M /M  in 
w e r e r 

terms  of  geometrical  variables.  The  question  arises  as  to  an  effective 
value  of  to  use  in  the  theory  for  non-circular  grains.  Since  the 
boundary  condition  is  being  applied  at  an  effective  radius,  a^,  a mass 
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balance  for  an  effective  radial  inflow  at  a^  yields 


2 n a t = p 
e 


V 

r 


S 

w 


or 


S 

w 

2tt  a ^ 
e 


(36) 


Equation  (35)  provides  a given  grain  and  Eq,  (36)  gives  a 

correction  to  to  be  used  in  the  theory  above,  Ihe  corrected 
physically  enters  to  guage  the  amount  of  flow  from  the  propellant  and 
its  effect  upon  the  acoustic  field.  For  large  area  grains  (e,g,,  a star 
compared  with  a circular  port)  lirg>  which  gives  the  same  physical 
effect  in  the  theory  as  occurs  in  practice;  that  is,  there  is  a stronger 
incoming  radial  flow.  Therefore,  the  approximations  of  this  section,  while 
somewhat  coarse,  are  expected  to  yield  proper  scaling  behavior  and  orders 
of  magnitude,  in  accord  with  the  purposes  here. 


Calculations  of  the  Green's  Function 

In  order  to  show  some  general  trends  and  to  see  where  Eq,  (34)  may 
be  a useful  approximation,  some  exact  calculations,  via,  Eqs,  (32)  and  (33), 
have  been  carried  out.  The  baseline  motor  parameters  used  in  computation 
are  shown  in  Fig,A-2,Mg=  0,3  is  a relatively  high  number  and  is  so  chosen 
in  anticipation  of  some  results  to  follow.  The  short  nozzle  assumption 
means  that  the  wavelength  under  consideration  is  long  compared  with  the 


nozzle  length.  Under  this  condition  the  nozzle  behaves  in  a quasi-steady 


fashion  and 


(20) 


M 

e 


It  may  readily  be  verified  in  the  theory  that  if  more  than  one  condensed 
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idiasc  group  is  present 


zd  = £ zd. 
i ^ 

where  is  calculated  on  the  basis  of  the  properties  of  the  i^^  group. 

It  is  assumed  for  these  calculations  that  A1  and  Al^O^  are  present  at 
an  average  mass  fraction  of  0,15  each  and  that  the  A1  diameter  is  300  p,m 
and  the  AI2O2  diameter  is  1 p,m.  Bits  yields  the  of^  and  02  values  in  Fig,  A^2, 

Two  curves  are  shown  in  Fig,  A-2,  One  is  an  exact  calculation  with  n*  0.5, 
and  one  is  the  approximate  solution  of  Eq,  (34),  There  is  generally  good 
agreement  in  trend  and  magnitude  except,  of  course,  the  approximate  solu- 
tion contains  no  resonance  phenomena.  The  error  in  the  approximate  solu- 
tion, however,  never  exceeds  10  dB  and  this  will  greatly  aid  in  order  of 
magnitude  arguments  later. 

The  resonances,  corresponding  to  the  first  and  second  longitudinal 
modes,  are  clearly  seen  in  Fig,  A-2,  The  nozzle  behaves  as  a relatively 
hard  wall,  and,  if  there  were  no  flow  present  the  resonance  of  the  first 
longitudinal  mode  would  occur  at  kl  = n*  The  flow  effect  shifts  the  actual 
resonance  frequency  slightly  down  from  kl  = tt. 

Figure  A-3  shows  the  Green's  function  for  various  pairs  of  values  of 
X and  y.  While  the  behavior  can  become  quite  complex  at  the  high  frequency 
end  of  the  spectrum,  below  about  200  Hz  the  relative  values  of  x and  y 
make  no  difference.  At  low  frequency  the  oscillation  gives  way  to  a nearly 
bulk  mode  oscillation  which  is  quite  well  predicted  by  the  approximate 
equation,  Eq.  (34),  It  should  be  noted  that  the  Green's  function  is  sym- 
metric; that  is  gy^(x,y)  = 

Figure  A-4  shows  the  effect  of  the  feedback  parameter^  for  various 
real  values  and  one  imaginary  value.  From  the  asymptotic  solution  of 
Eq,  (34)  it  is  seen  that  jT  = (5*^)/?  is  a critical  value  for  the  low 
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Figure  A-3.  Green's  Function  for  Various  Pairs  of  Values 
of  X and  y. 
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frequency  solution..  By  calculation,  at  the  resonance  point  of  the  first 
longitudinal  mode  'JJ's  2,36  is  a critical  value.  Above  these  values,  and  if 
iT  is  real,  a stationary  random  oscillation  is  impossible  and  an  unstable 
oscillation  will  occur.  It  is  assumed  here  that  jT  is  always  below  the 
critical  stability  value.  The  imaginary  value  of~  primarily  effects 
the  resonant  frequency  at  moderate  frequencies,  but  at  low  frequency 
it  also  affects  the  magnitude  of  g^,  as  may  be  seen  from  Eq,  (34)  and 
Fig,  A-4, 

InFig,A-5  the  effect  of  Mach  number  is  shown.  Both  M and  must  be 
varied  simultaneously  for  a given  motor  geometry  as  demanded  by  Eq,  (35), 
The  Mach  number  is  varied  by  a factor  of  two  in  Fig,  A-5,  The  effect  of 
lowering  M is  to  decrease  the  nozzle  damping.  This  is  especially  effective 
at  low  frequency  and  at  resonance  points. 

In  Fig,  A»6  the  effect  of  M^  alone  is  shown  under  the  constraint  that 
the  motor  mass  flow  (thrust)  is  the  same.  By  Eq,  (36),  this  constraint 
requires  M^/a  to  be  maintained  constant.  Consequently,  Fig,  A-6  is  inter- 
preted as  the  effect  of  a change  in  grain  configuration  to  accommodate  a 
slower  burning  propellant.  There  is  an  important  effect  of  this  variable. 
Slower  burning  propellants  raise  g^. 

Finally,  the  effect  of  particle  damping  is  shown  in  Fig,  A-7,  There  is 

a surprisingly  small  effect  when  the  particle  mass  fraction  is  changed, 

(22) 

It  was  anticipated  by  example  in  the  Instability  field  that  the  effect 
would  be  stronger.  In  fact,  here  the  effect  is  negligible.  If  one  reviews 
the  equations  it  can  be  seen  that  the  comparative  magnitude  of  z^  terms 
is  quite  small. 

It  appears  that  only  near  unstable  resonance  points  can  the  effect 
show  up  strongly.  This  is  limited  by  the  behavior  at  the  longitudinal 
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Figure  A-5.  Green's  Function  for  Two  Values  of  Exit 
Mach  Number. 
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mode  resonance  point.  By  and  large,  however,  the  particle  damping  effect 
is  negligible  for  these  stationary  random  oscillations  for  moderate  jj,. 
Not  shown  are  the  effects  of  particle  size  change,  because  they  are  as 
ineffective  as  mass  fraction  changes  in  affecting  the  curves, 

summarizing,  the  behavior  of  the  g^  function  is  strongly  affected 
by  feedback  parameter,  flow  speeds,  and  bum  rate,  but  is  insensitive 
to  particle  effects.  The  low  frequency  solution  is  adequate  for  order 
of  magnitude  arguments  sufficiently  far  away  from  onset  of  the  first 
longitudinal  mode. 


Decomposition  of  the  Solution 

Viewing  Eq,  (28),  it  is  desired  to  further  manipulate  the  solution 
into  a form  suitable  for  computation.  Consider  first  the  first  term  on 
the  right  hand  side  of  Eq,  (14),  Repeated  application  of  the  divergence 
theorem,  using  the  assumption  that  = 0 on  all  S except  and  using 
Eq,  (26),  the  following  operations  are  carried  out,: 
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1 5 
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I { ty,  L®a  ly  ] ' 

f7  J i 


^ o J i;  J ^ 
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c 2 
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S ^ S 

e e 


dV 


111 


' J L®«.  I?,  ] ■‘S  + J [ + (2  + *d)  ‘•“el  ■*= 

s ^ s 

e e 

- •'iiO'i)  - x'  1 

V (37) 

In  developing  Eq.  (37),  terms  of  the  order  of  compared  to  unity  have 
been  dropped,  consistent  with  the  procedure  used  for  the  right  hand  side 
of  Eq,  (l4).  The  first  term  on  the  right  hand  side  of  Eq,  (37)  is  a 
rather  complicated  interaction  of  the  turbulence  with  the  nozzle  entrance 
plane.  The  next  three  terms  only  depend  upon  with  the  first  being 
a nozzle  entrance  plane  interaction  the  second  being  the  locally  observed 
density  fluctuation  due  to  the  turbulence  and  the  third  is  a volume  dis- 
tribution effect.  Similar  manipulation  can  be  made  Involving  the  term 
which  comprises  F,  The  second  term  of  the  right  l.ind  side  of  Eq,  (28)  is 
obtained  by  plugging  in  the  appropriate  p values  on  the  appropriate  sur- 
faces, There  is  then  a reasonable  amount  of  cancellation  which  occurs 
among  various  terms.  The  final  solution  for  appears  as 

^(pi'abcab  ab 
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II  = - ikp  M 1 ???  g dS 

a r J u)  '’uj 
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“b  ' - f Pe  J ^ i ■*'' 

® V 

III  = ikp  f c b g dS 

J UJ  u)  u) 

S 

e 

IV  = - ^ (T.,) 

a c 11  (11 

e 

IV,  = - D O /^ox 

b Pe  (38) 

Notice  also  Eq,  (13)  so  that  if  one  calculates  the  pressure  transform, 

term  IV does  not  enter  the  pressure  calculation.  Moreover,  at  the  head 

end  and  side  walls  IV  =0  so  that  terms  I -III  are  the  only  terms  of 

a 

interest  in  the  wall  pressure,  and,  hence  the  vibration  problem. 

Discussion  of  Eqs.  (38)  is  now  in  order.  Term  I is  an  effect  of 

A 

distributed  turbulence  on  pumping  up  axial  oscillations.  It  depends  only 

on  the  axial  components  of  Tjj,  Term  Ij^  is  an  interaction  of  the  axial 

Tji  with  the  nozzle.  So  is  term  I , Terms  II  and  II.  are  direct  com- 

cab 

bustion  noise  terms,  with  II^  caused  near  the  propellant  surface  and  IIj^ 
by  burning  of  any  particulate  matter  in  the  bulk  volume  of  the  chamber. 

Term  III  is  an  interaction  of  hot  spots  with  the  nozzle  and  is  called  en- 
tropy noise.  Note  it  is  directly  proportional  to  the  entropy  admittance 

coefficient,  b . 

tu 

This  solution  is  actually  quite  remarkable  in  its  simplicity  and  cap- 
ture of  the  essential  physics.  It  contains  all  effects  which  were  antici- 
pated at  the  outset.  This  was  not  guaranteed  because  the  approach  is  an 
approximate  one  and  the  Lighthill  approach  has  not  been  previously  attempted, 
to  the  author's  knowledge,  on  an  interior  acoustics  problem.  In  independent 
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calculations,  not  presented  here,  the  term  which  will  survive  of  terms 

and  turns  out  to  be  identical  with  a vorticity-nozzle  interaction 

term  derived  if  the  formulation  used  is  the  linearized  acoustic-vorticity 

approach  to  aeroacoustics.  Term  111  is  similar  in  form  and  scales  in  the 

(12) 

same  manner  as  entropy  noise  derived  from  a different  approach.  Term 
11^  is  identical  in  form  with  previous  results  in  combustion  noise^^  ^ 
Therefore,  there  is  a high  degree  of  confidaice  that  the  solution  of 
Eqs,  (38)  represents  a complete  picture  of  the  actual  solution,  will 
yield  proper  orders  of  magnitude  and  will  yield  proper  scaling  laws. 

The  problem  now  is  to  obtain  numbers  to  set  into  the  theory  to  make  the 
actual  computations. 
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